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Abstract: We study the singular values of the Dirac operator in dense QCD-like theories 
at zero temperature. The Dirac singular values are real and nonnegative at any nonzero 
quark density. The scale of their spectrum is set by the diquark condensate, in contrast to 
the complex Dirac eigenvalues whose scale is set by the chiral condensate at low density and 
by the BCS gap at high density. We identify three different low-energy effective theories 
with diquark sources applicable at low, intermediate, and high density, together with their 
overlapping domains of validity. We derive a number of exact formulas for the Dirac singular 
values, including Banks-Casher-type relations for the diquark condensate, Smilga-Stern- 
type relations for the slope of the singular value density, and Leutwyler-Smilga-type sum 
rules for the inverse singular values. We construct random matrix theories and determine 
the form of the microscopic spectral correlation functions of the singular values for all 
nonzero quark densities. We also derive a rigorous index theorem for non-Hermitian Dirac 
operators. Our results can in principle be tested in lattice simulations. 
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1 Introduction 



A prominent feature of quantum chromodynamics (QCD) in the vacuum is the dynamical 
breaking of chiral symmetry through the formation of a chiral condensate (ipip) ■ While an 
analytical demonstration of this phenomenon from the underlying theory is still lacking, 
a number of studies have deepened our understanding of its nature. For example, spon- 
taneous chiral symmetry breaking is reflected in the infrared limit of the Dirac eigenvalue 
spectrum through the Banks-Casher relation [1] and the Smilga-Stern relation [2]. Also, a 
universal finite-volume domain (called the e-regime) has been discovered [3] in which the 
theory becomes zero-dimensional and is governed entirely by the global symmetries of the 
system. It was shown to lead to an infinite number of constraints on the low-lying Dirac 
eigenvalues, known as the Leutwyler-Smilga spectral sum rules [4]. 

Thanks to the universality, rigorous results for the spectral properties of the Dirac op- 
erator can be derived from a much simpler chiral random matrix theory (RMT) which has 
the same global symmetry-breaking pattern as the QCD vacuum (see [5, 6] for reviews). 
This has advanced our knowledge of the low-lying Dirac eigenvalues to the level of micro- 
scopic spectral correlation functions, which eventually made it possible to determine the 
value of the chiral condensate to high precision by first-principle lattice QCD simulations 
[7], thus numerically verifying that chiral symmetry is spontaneously broken in the QCD 
vacuum. There are also other QCD-like theories that can be described by RMT. They can 
be classified by the Dyson index f3. If the Dirac operator commutes with an anti- unitary 
operator T, then f3 = 1 if T 2 = 1 and f3 = 4 if T 2 = -1. Otherwise /3 = 2. Equivalently, 
(3 = 1, 2, or 4 if the representation of the gauge group in which the fermions transform is 
pseudoreal, complex, or real, respectively. 1 

At large quark chemical potential /i, 2 QCD is believed to exhibit another nonper- 
turbative phenomenon: The attractive interaction between quarks near the Fermi surface 
leads to color superconductivity characterized by a diquark condensate {i^ip) (see [9, 10] 
for recent reviews). In the color-superconducting phases, gluons acquire masses through 
the Anderson-Higgs mechanism, and weakly interacting quarks acquire a Bardeen-Cooper- 
Schrieffer (BCS) gap A. In particular, in the color-flavor-locked (CFL) phase [11] of Nf = 3 
QCD at high density, chiral symmetry is spontaneously broken in much the same way as 
in the QCD vacuum, but the order parameter is the four-quark condensate rather than the 
conventional chiral condensate. Using effective-theory techniques, it was shown [12] that 
the distribution of the smallest Dirac eigenvalues in the CFL phase is entirely governed by 
global symmetries and that the relevant scale of the low-lying spectrum is set by the BCS 
gap A. 

Unfortunately, at /i ^ the complex phase of the fermion determinant in the QCD 
partition function has so far hampered lattice computations of physical quantities such as 
the BCS gap and the diquark condensate. This is an example of the so-called sign problem, 

1 This can be shown following the considerations in sections 21.2-21.4 of [8]. 

2 Throughout most of the paper we will use the terms "chemical potential" and "density" interchangeably, 
unless the distinction is important. Nonzero density always implies nonzero /x. However, the reverse is not 
always true, e.g., in two-color QCD at small the density is zero for fi < m^/2, where m n is the pion mass. 
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which is encountered in many areas of physics. However, several QCD-like theories are 
devoid of the sign problem even at nonzero /i [13] and develop a nonzero BCS gap and 
a diquark condensate (or pion condensate) signaling superfluidity. Examples include two- 
color QCD (/3 = 1), Sp(2iV c ) gauge theories with fundamental fermions and an arbitrary 
number N c of colors (/3 = 1), SU(iV c ) gauge theories with adjoint fermions (/? = 4), and 
SO(N c ) gauge theories with fundamental fermions (/3 = 4). Another example is three- 
color QCD at nonzero isospin chemical potential ((3 = 2). 3 In spite of the clear differences 
between these theories and three-color QCD at nonzero quark chemical potential, they 
share some common features with bona fide QCD, e.g., the same mechanism of quark- 
quark pairing that leads to superfluidity or color superconductivity at nonzero ^, and the 
universality of their phase diagrams [14]. This observation leads us to expect that we 
can obtain some insights into the properties of realistic QCD matter through a deeper 
understanding of these QCD-like theories. 

Indeed, the aforementioned work on the Dirac eigenvalues in the CFL phase has been 
successfully generalized to two-color QCD at high density (/3 = 1) [15], where the BCS pair- 
ing was shown to have sizable consequences on the behavior of the small Dirac eigenvalues, 
as represented through new Leutwyler-Smilga-type sum rules. (This analysis permits a 
straightforward extension to theories with (5 = 2 and 4.) Moreover, the chiral random 
matrix theory that describes the Dirac eigenvalue distribution on the scale characterized 
by the BCS gap was constructed [16] and solved analytically [17, 18]. These results make 
it possible in principle to compute the BCS gap from Dirac spectra on the lattice, which 
would directly verify the BCS-type superfluid phase of QCD-like theories. 

Although the BCS gap A and the diquark condensate (tpip) are closely intertwined 
at high density, they are a priori different objects. While it is unclear to what lower 
densities the BCS-type pairing with well-defined A will persist, the diquark condensate 
is still nonvanishing and breaks global symmetries even at lower densities in those QCD- 
like theories. Some of the Nambu-Goldstone (NG) modes are diquarks charged under 
baryon number (with mass m^) so that they exhibit Bose-Einstein condensation (BEC) 
for n > [14, 19-25]. 4 The Bose-Einstein condensate at small /U consists of strongly 

coupled bound diquarks, in contrast to the weakly coupled BCS-type diquarks at large /i. 
Despite this qualitative difference, however, the quantum numbers of the condensate are 
the same at small and large //, so it is natural to expect that there is no phase transition 
between the two limits. This phenomenon is called (relativistic) BEC-BCS crossover. It 
passes a number of nontrivial tests [21, 22, 26], and its possible realization in dense quark 
matter in QCD has been investigated in various model calculations [27-42] (see also [43-48] 
for related models). In view of the known relation between A and the Dirac eigenvalues 
[12, 15], a natural question arises: How is the diquark condensate {ipvp} reflected in the 

3 To evade the sign problem, one has to assume an even number of flavors with degenerate quark masses 
in the ft = 1 and /3 = 2 cases. On the other hand, we can take an arbitrary number of flavors with 
nondegenerate quark masses in the /3 = 4 cases. 

4 In QCD at nonzero isospin chemical potential fii, the usual pions, which do not have baryon charge 
but isospin charge, show BEC, i.e., (d'y^u) ^ 0, for \ii > m^. In this case, however, it is still possible to 
view this as a BEC of "diquarks" if we switch d — >• d' = Cd . 
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Dirac spectrum of QCD at nonzero density? 

In this paper, we point out that it is the spectrum of the singular values of the Dirac 
operator D(fi), i.e., the square roots of the eigenvalues of D(fj,y D(/j,), that carries the 
information on the diquark condensate (i/j^) at any nonzero chemical potential in all QCD- 
like theories we consider here. By inserting a diquark source (instead of a quark mass) 5 we 
derive a number of rigorous relationships covering the entire BEC-BCS crossover region, 
such as Banks-Casher-type relations, Smilga-Stern-type relations, and Leutwyler-Smilga- 
type spectral sum rules for the Dirac singular values. This significantly extends previous 
related work [49] in many directions. We also construct chiral random matrix theories and 
determine the form of the microscopic spectral correlation functions for singular values at 
all nonzero densities. In most of this paper we will work in the chiral limit. 

Note that, although the Dirac singular values and the Dirac eigenvalues coincide at 
fi = 0, they are essentially different objects at nonzero /U. 6 First, while the Dirac eigenvalues 
are complex, the Dirac singular values are always real and nonnegative. Second, while the 
scale of the small Dirac eigenvalues is governed by the chiral condensate at small \i and 
by the BCS gap at large /i, the scale of the small Dirac singular values is governed by 
the diquark condensate at any \x. From the viewpoint of symmetries, the diquark source 
breaks U(l)# (baryon number), and therefore the spectrum of the Dirac singular values 
characterizes the U(l)# symmetry breaking. This is in contrast to the spectrum of the 
Dirac eigenvalues, which is unrelated to U(l)s since neither the quark mass (the source 
for the chiral condensate at low density) nor the quark mass squared (the source for A 2 at 
high density) breaks U(l)s- 

Our results for the Dirac singular values will make it possible to measure the magnitude 
of the diquark condensate with high precision by lattice QCD simulations at any density. 
Together with earlier results for the Dirac eigenvalues they lead to a more detailed under- 
standing of superfluidity and the BEC-BCS crossover in QCD-like theories, and hopefully 
also of the physics of color superconductivity in three-color QCD. 

The structure of this paper is as follows. In section 2, we introduce the microscopic 
theories considered in this paper and write down the general expressions for the partition 
functions. In section 3, we review basic properties of the eigenvalues and singular values of 
the Dirac operator and clarify their relation at nonzero chiral chemical potential. Special 
emphasis is given to the zero modes. In section 4, we derive Banks-Casher-type relations 
for two-color QCD (/3 = 1), QCD at nonzero isospin chemical potential (/3 = 2), and QCD 
with adjoint fermions (/3 = 4). In sections 5, 6, and 7, we concentrate on two-color QCD 
(J3 = 1) for simplicity and brevity. The results of these sections admit straightforward 
generalization to the other QCD-like theories with (3 = 2 and 4. In section 5 we introduce 
three different low-energy effective Lagrangians with diquark sources applicable at low, in- 
termediate, and high density. In sections 6 and 7, we derive Smilga-Stern-type relations and 

5 The diquark source can also be called Majorana mass while the quark mass is called Dirac mass. 

6 As an aside we mention that at fj, — the spectrum of the so-called Hermitian Wilson-Dirac operator in 
lattice QCD with Wilson fermions is nothing but the singular value spectrum of the Wilson-Dirac operator 
owing to the 75-Hermiticity of the operator. In this case, the density of the singular values near the origin 
is proportional to the parity-breaking condensate [50]. 
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Leutwyler-Smilga-type sum rules, respectively, using the effective theories from section 5. 
In section 7 we also construct chiral random matrix theories that describe the spectrum 
of the Dirac singular values in the e-regime and determine the form of the microscopic 
spectral correlation functions. Section 8 contains the conclusions and an outlook. 

In appendix A, we summarize our definitions and conventions. Appendix B describes 
the derivations of the singular value representations of the partition functions for the theo- 
ries with P = 1, 2, and 4. In appendix C, we comment on the importance of the positivity 
of the measure, which could be spoiled by the diquark sources. We also discuss QCD 
inequalities and derive constraints on the symmetry-breaking pattern for positive definite 
measure. In appendix D, we present careful derivations of the anomaly equation and the 
extension of the index theorem to /U ^ 0, paying special attention to the non-Hermitian 
nature of the Dirac operator. In appendix E the derivation of eq. (6.19) in the main text 
is outlined. In appendix F we comment on the random matrix theory for QCD at nonzero 
isospin chemical potential (/3 = 2). Finally, appendix G is devoted to the derivation of 
eq. (7.52) in the main text. 

2 Microscopic theories 

In this section, we introduce QCD-like theories with the Dyson indices j3 = 1, 2, and 4, 
emphasizing the anti-unitary symmetries of the Dirac operator and the global symmetries 
of the theories. Since the analysis of all three cases is similar, we first examine the f3 = 1 
case in detail and then discuss /3 = 2 and /3 = 4 briefly without redundancy. We take 
two-color QCD and QCD with adjoint fermions as examples for the /3 = 1 and /3 = 4 cases, 
respectively. The same arguments are readily applicable to Sp(2iV c ) gauge theory (/? = 1) 
and SO(N c ) gauge theory (/3 = 4). 

Unless stated otherwise we always work in Euclidean space and at zero temperature. 
Our definitions and conventions are summarized in appendix A. 

2.1 Two-color QCD {(3 = 1) 

The Dirac operator of two-color QCD in the presence of a chemical potential is given by 

D(ji) = ^ U D U + ^74 with D u = d v + %A V , (2.1) 

where A v = A^T a /2 is the gauge field and the r a are the generators of SU(2) color, i.e., the 
Pauli matrices. For fi = the Dirac operator is anti-Hermitian, but for \x ^ it no longer 
has definite Hermiticity properties since /i74 is Hermitian. However, for two colors there is 
an anti-unitary symmetry that can be expressed in two equivalent ways [4, 19, 51, 52], 

[Ct 2 K, iD(fi)} = with {Ct 2 K) 2 = t , (2.2a) 
[j 5 Ct 2 K, D(h)} = with {^Ct 2 K) 2 = 1 , (2.2b) 

where C = i^Al2 is the charge conjugation matrix (see appendix A) and K is the operator 
of complex conjugation. Hence the Dyson index is (3 = 1 in this case and we can choose 
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a basis in which the Dirac operator is real. This implies that det D(/i) is real, which also 
follows from Ct275-D(//)Ct275 = D(fi)* as a consequence of (2.2b). 

We introduce Nf quark flavors 7 described by Dirac spinors ipf (/ = 1, . . . , Nf), which 
we collect in tp = (ifii, . . . ,ipN f ) T . Each Dirac spinor can be split into two Weyl spinors, 
which we collect in tpR = (ipm, . . . ,ipN f R) T and = (ipiL, ■ ■ ■ ,^N f L) T - The fermionic 
part of the Lagrangian, including mass term and diquark sources, is given by 

C f = $[D(ji) + MP L + M^P R ]iP + ^ t Ct 2 (JrPr + JlPl)^ + \^Ct 2 {J r Pr + j\P L W, 

(2.3) 

where Pr/l = (1 ±75)/2 and M, Jr, and Jl are A^y-dimensional complex matrices in 
flavor space. As a consequence of the Pauli principle, it suffices to take Jr and Jl to 
be antisymmetric since their symmetric parts drop out of the combinations in (2.3). For 
greater generality we have allowed for two independent matrices Jr and Jl. 

We briefly summarize the symmetries of (2.3) obtained in [19, 20], assuming M = ml 
and Jr = — Jl = jl, where m and j are real and I is defined in (A. 9). For fj, = 0, the 
Lagrangian in the absence of mass term and diquark source is symmetric under SU(2./Vy). 8 
The mass term breaks this symmetry to Sp(2iVy). The diquark source transforms into the 
mass term under an SU(2iVy) rotation, and therefore it brings nothing new at fi = 0. For 
/x ^ and no sources, the chemical potential breaks the SU(2iVj) symmetry explicitly to 
SV(Nf) R x SU (N f ) l x U(l) b- The mass term breaks this symmetry to S\J(N f ) L+R x U(1) B , 
while the diquark source breaks it to Sp(iV/)# x Sp(iV/)i. In the presence of both mass 
term and diquark source the symmetry is broken to Sp(Nj)l+r- 

We now express Cf in the so-called Nambu-Gor'kov formalism. Defining 




* = _r T , (2.4) 



(2.3) can be rewritten as 



C f = ^ T W^> (2.5) 



with 



w = | Ct 2 (J r Pr + J L P L ) -D(v) t - M T P L - M*P R \ , , } 

1 D(p) + MP L + M^Pr -Ct 2 (J r P l + J ] l Pr) 1 ' 



Since Jr, Jl, C, and t 2 are all antisymmetric, W is antisymmetric as well, and thus 

J exp \-\J d * x ^ T W^ = Pf (W) , (2.7) 
7 We assume Nf < 11 to ensure asymptotic freedom. 

8 U(l)s is contained in SU(2iV/). There is an additional U(1)a symmetry which is anomalous and 
therefore not considered here, but see section 5.3. 
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where Pf denotes the Pfaffian, which for an antisymmetric matrix of even dimension N is 
defined as 9 

Pf (X) = 2^/ 2 (N/2)l ^ s S n ( (J ) X -(i)-(2) • • • x *(N-i)*(N) ■ (2.8) 
The partition function is therefore 

Z =( P W> YM ' ^ 

where the subscript YM means that the average is over gauge fields weighted by the pure 
gauge (Yang-Mills) action. Since in this paper we will almost always work in the chiral 
limit we set M to zero in (2.6). As shown in appendix B.l, we then obtain from (2.9) 



z(j L , j r ) = ( [ pf ( j r ) pf ( 4)] nR [ pf (4) pf ( j l ) 



det'(D(fi)W(fi) + J R J R P R + jlj L P L )) , (2.10) 

/ YM 

where n R (til) denotes the number of right- (left-) handed zero modes of D^D and the prime 
on the determinant means that the zero modes of D^D are omitted. 10 This expression is 
invariant under S\J(Nj) R x S\J(Nf)i, as it should be, since Pf(J R ) — > Pf(Uj,J R U R ) = 
det(U R ) Pf (Jr) = Pf(J fl ) and likewise for Pf(J L ). 

We can also add a P- and CP- violating term i9F F /32n 2 to the Lagrangian, where 
a/3 = ^ £ a/3^8^s- ^he topological charge v of a gauge field configuration is given by 

VsiftA (2.11) 



32vr 2 



and related to the number of zero modes of D by v = n R — n^. 11 The #-term corresponds, 
for a fixed gauge field, to a term ivQ in the action. The partition function is then given by 

oo 

Z{9)= e "^> ( 2 - 12 ) 

v=— oo 

where Z v is computed by integrating only over gauge fields with fixed topology v. 

For zero diquark sources and nonzero quark masses it is well known that an axial 
rotation of the quark fields changes the fermionic measure in such a way that the #-term 
can be traded for a redefinition of the quark mass matrix, M M e~ ie / N f [4]. 12 For 
nonzero diquark sources and in the chiral limit, it follows from (2.10) that the #-term can 



9 Later we will use the following basic properties of the Pfaffian: Pf(X) 2 = det(X), Pf(X T ) = 
(-l) N/2 Pf(X), Pf(cX) = c N/2 Pf(X) for c G C, and Pi(UXU T ) = det((7) Pf(X). 

10 The numbers ur and ul are also equal to the number of right- and left-handed zero modes of D, see 
the discussion in section 3. 

11 Note that for (i / the equality v — ur — ul is violated on a gauge field set of measure zero, see 
appendix D. Here (and also in the next two subsections) we exclude this possibility. 

12 In [4] M — > M e l6 / N i is used, which is due to different conventions, see also footnotes 52 and 69. 
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be traded for a redefinition of the diquark sources, Jr — > Jr e t6 Rl N f and Jl — > Jl^ Ql ^ s 
with 9 = {Or - 9 L )/2. 

Note that for some choices of the diquark sources and/or for a nonzero value of 6 
the fermionic measure in the partition function is not positive definite, which causes some 
subtleties that are discussed in detail in appendix C. 

2.2 QCD with isospin chemical potential {(3 = 2) 

We now consider Nf = 2 QCD at nonzero isospin chemical potential pi = 2p for an 
arbitrary number of colors N c > 3. The Dirac operator is given as in (2.1), except that 
the r a are now the generators of SU(iV c ) in the fundamental representation. The chemical 
potential p is assigned to the u-quark, while —p is assigned to the (i-quark. For p ^ 0, the 
Dirac operator is no longer anti-Hermitian, but because of D(py = —D(—p) the fermion 
determinant is real and nonneg ative: det D(p)D(-p) = det D(ji)D{jny* > 0. 13 Since there 
is no anti-unitary symmetry we have (3 = 2. 

The fermionic part of the Lagrangian is given by 

£ = u(~f u D u + p^)u + d(^ u D u - /i^4)d + (mu R UL + md R d,L + h.c.) 
+ (X'u^cLr + pd\un + h.c.) 

/_ j\ ( D(ji) + mP L + m*P R \*Pr + p* P L 
V )\ pP R + XPl D{-p) + mP L + m*P R 

where m is the degenerate mass of the two quarks, and p and A are "pionic" sources. From 
this we obtain (see appendix B.2) 

Z(p, A) = ((- P A*) nR (-p*\) nL det'(rtD + pp*P R + AA*P L )) yM , (2.14) 

where the prime again indicates that the zero modes of D^D are omitted. A nonzero 
#-term can be introduced in (2.14) by redefining the pionic sources as p — > pe l9R ^ 2 and 
A -»• Xe i0L / 2 , where again 9 = (Or- 6 l )/2. 

The symmetries of (2.13) with p = —A € 1R are as follows. The Lagrangian at p = 
in the absence of sources is symmetric under SU(2)i x SU(2)# x U(1)b, which is broken to 
SU(2) x U(1)b by the mass term or the pionic sources. 14 There is also a U(1)a symmetry 
that is broken by the anomaly. Without sources the chemical potential breaks the SU(2)x, x 
SU(2)ij x U(1)b symmetry to U(1)l x XJ(1)r x U(1)_b, where 11(1)^^ is generated by 
ts G su(2) L / R . The mass term breaks this symmetry to U(1)l+r x U(1)b, while the 
pionic sources break it to \J(1)l-r x U(1)_b- With both mass term and pionic sources the 
remaining symmetry is U(l)#. 

13 Note that det(— D^) = det D* if has no zero modes. If it has zero modes we simply have = 0. 
14 For the mass term we end up with SU(2)_l+_r, while for the pionic sources we end up with a different 
SU(2) subgroup given by the condition \J\XiUl = ti, where ti is the second generator of SU(2). 
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2.3 QCD with adjoint fermions (/3 = 4) 

Finally we consider QCD with fermions in the adjoint representation. 15 The Dirac operator 
in the presence of a chemical potential is given by 

D(ji)= 7„A, + /*74 with D v = d u 8 ab + (f) ab A c u , (2.15) 

where (f c ) a b = fabc denotes the generators of the adjoint representation (or structure 
constants). For fi = the Dirac operator is anti-Hermitian, but for [i ^ it again loses its 
Hermiticity properties. There is an anti-unitary symmetry that can be expressed in two 
equivalent ways [4, 19, 51, 52], 

[CK,iD(fi)] = with {CKf = -1 , (2.16a) 
[l 5 CK,D(v)] = Q with ( 75 CK) 2 = -1. (2.16b) 

Hence /3 = 4 and we can choose a basis in which the Dirac operator is quaternion real. The 
symmetry (2.16b) implies C"f^D{^i)C^^ = D(/i)*, from which it follows that det D(fi) is real 
(and actually nonnegative because all eigenvalues occur in quadruplets, see section 3.1). 
Because of (CK) 2 = — 1 one can show that for \x = (but not for \i ^ 0) the eigenvalues 
of the Dirac operator are twofold degenerate with linearly independent eigenstates ip and 
Cip* (Kramers degeneracy). 

Because the adjoint representation is real, it may be convenient to describe the fermions 
in the partition function in terms of Majorana fields. However, Majorana fermions cannot 
be defined in Euclidean space, and therefore we first write the Lagrangian in Minkowski 
space and then analytically continue to Euclidean space by a Wick rotation [4, 53]. In 
Minkowski space, the Lagrangian for Nf = 1 Dirac fermions with diquark sources reads 

4/ /=1) = $(rfDv ~m- ^7 V + \ WC{j R P R + j L P L )V + h.c] 

-C(j* R P L + 3lPr) il v D v - m - M7° 
C(-iYD v + m - ^°)C CUrPr + jlPl) 

C(j* R p L + 3l P R) ~C{il v D v -m- ^7°) 
C(-iYD u +m- ^7°) C(j R P R + j L P L ) 





(2.17) 



— T 

where in the last line we have defined ip c = Ctp . The partition function is thus given by 

7 (iV/=i) _LJ C(f R P L + j* L P R ) -Cii-fD v - m - M7° A \ (l? 1 s . 

" \ \C(-vfD v + m - /x 7 °) C(j R P R + 3L P L ) )/ ym (2 - 18j 

This result generalizes trivially to general Nf. The partition function is then simply the 
expectation value of a product of Nf Pfaffians, possibly with different masses. 

For simplicity we now take the chiral limit. After Wick rotation, the partition function 
for Nf flavors in Euclidean space becomes 



c( j* r p l + r L p R ) CD(fx) 

-CD(^ C(J R P R + J L P L ) i 



z E = ( pf ( y R J:Z .X H> „, r :T, - I > . (2-i9) 



YM 



3 The results obtained for the adjoint representation easily generalize to other real representations. 
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where Jl and J R are now symmetric matrices of dimension Nf. For (5 = 4, Nf does not 
have to be even. In appendix B.3 we show that this leads to 

Z(J L , Jr) = (det(- J R j{) nR/2 det(-J R J L ) n ^ 2 det"(D j 'D + J R J R P R + J ] l J l P l )) ym , 

(2.20) 

where det" indicates that the zero modes of D^D are omitted and that each degenerate 
eigenvalue of D^D is counted only once. It follows from the derivation in appendix B.3 
that a nonzero #-term can be introduced in (2.20) by redefining J R — > J R e %eR / 2N f Nc and 
Jl -> J L e WL/2N f Nc , where again 6 = (6 R - 6 L )/2. 

The symmetries of (2.17) extended to general Nf with degenerate masses and J R = 
—Jl = jl with real j are as follows [20]. The Lagrangian at /x = in the absence of 
sources is symmetric under SXJ(2Nf), which is broken to SO(2Nf) by the quark mass or 
the diquark source. 16 There is also the usual anomalous \J(1)a symmetry. Without sources 
the chemical potential breaks the SU(2iV/) symmetry to SU(JV/)i x S\J(Nf) R x U(l)s. 
The mass term breaks this symmetry to SU (Nf)i +R x U(1)b, while the diquark source 
breaks it to SO(Nf)i x SO(Nf) R . With both mass term and diquark source the remaining 
symmetry is SO(N f ) L+R . 

3 Eigenvalues and singular values of the Dirac operator 

In this section we discuss the eigenvalues and singular values of the Dirac operator and 
related quantities. Two preliminary remarks are in order. 

First, some parts of the discussion rely on the index theorem. In appendix D we show 
that for a non-Hermitian Dirac operator such as -D(/i) the index theorem takes the form 

^ J d A xFF = ^[mdD( f i) + mdD(^] (3.1) 

and that 

ind D(fi) = ind D{ny almost surely, (3-2) 

where mdD(fi) = dimkerD R — dimkerD^, see (D.l) for the notation. The meaning of 
"almost surely" is that to have indD(/i) ^ indD(/i)' the gauge field needs to be fine-tuned, 
which corresponds to a gauge field set of measure zero. In this section we ignore this set 
of measure zero and use the index theorem in its usual form. 

Second, we implicitly assume a regularization (such as lattice QCD) that allows us to 
count the number of eigenstates. Some of our arguments rely on the index theorem, which 
can be (and usually is) violated by the regulator. Therefore our results apply only after 
the regulator has been removed. We assume that the procedure of removing the regulator 
does not invalidate the results that rely on the index theorem. 

16 As in the case of /3 = 1, the diquark source transforms into the mass term under an SU(2JV/ ) rotation. 
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3.1 Eigenvalues 

The discussion of this section extends that of [52]. The eigenvalue equation for the Dirac 
operator is 

D(fi)lp n = \ n ip n . (3.3) 

For simplicity we will omit the argument \x if no confusion is likely to arise. 

Let us first consider /j = 0. 17 In that case the eigenvalues A n are purely imaginary, 
and because of {.0,75} = the nonzero eigenvalues come in pairs ±A n . If the eigenstate 
corresponding to A n is ip n , then the eigenstate corresponding to — A n is 7sVVi- There can 
also be eigenvalues equal to zero, and the corresponding eigenstates can be arranged to be 
eigenstates of 75. In general there are tir (ul) right-handed (left-handed) zero modes with 
eigenvalue +1 (—1) of 75. The difference tir — til is equal to the topological charge v of the 
underlying gauge field configuration via the index theorem and stable under perturbations 
of the gauge field. All other possible zero modes are accidental in the sense that they require 
a fine-tuning of the gauge field. This implies that generically there are only zero modes of 
one chirality. In the remainder of this section we assume that there are no accidental zero 
modes. 

We now discuss what happens to the eigenvalues as fi is turned on (for a fixed gauge 
field and in a finite volume) and first consider (3 = 1 (e.g., two-color QCD). Since D(fi 7^ 0) 
is no longer anti-Hermitian one would expect the eigenvalues to move into the complex 
plane for any [i. However, using the symmetry (2.2b) one can show that the nonzero eigen- 
values come either in quadruplets A, — A, A*, — A* (with eigenstates ifi, 75^, Ct2^i(>* , Ct2?/>*) 
if A is complex or in pairs ±A (with eigenstates tp,^^i(j) if A is purely real or purely imag- 
inary. Since at fj, = the nonzero eigenvalues are generically nondegenerate (because of 
level repulsion due to interactions) , a quadruplet cannot be formed for infinitesimally small 
/x. 18 What happens (see figure 1) is that as a function of /u, the eigenvalues move along the 
imaginary axis until two eigenvalues (and their partners) become degenerate, i.e., the effect 
of [i overcomes the level repulsion. As /x is increased further, these four eigenvalues move 
into the complex plane and form a quadruplet. Another possibility is that a pair of origi- 
nally imaginary eigenvalues hits zero and then becomes a pair of real eigenvalues. Finally, 
two real eigenvalues (and their partners) can also merge and then become a quadruplet. 19 
The original topological eigenvalues equal to zero (but not the accidental ones) are stable 
under the perturbation by /i. The corresponding zero modes change smoothly with [i and 
remain eigenstates of 75 [54-57]. 

For (3 = 2 there is no anti-unitary symmetry, and the eigenvalues (which still come in 
pairs ±A with eigenstates ift, 75 V>) move into the complex plane immediately as \i is switched 
on. Exactly real or imaginary eigenvalues only occur accidentally, i.e., if the gauge field is 
fine-tuned. The topological zero modes behave as in the case of (3 = 1. 

17 The arguments for fi = apply to all values of the Dyson index f3, except that for ft — 4 we have the 
additional Kramers degeneracy already mentioned in section 2.3. 
18 The assumption of a finite volume is essential here. 

19 A11 of these three possibilities have been verified in lattice simulations and random matrix studies. We 
thank Jacques Bloch for performing the lattice simulations. 
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Figure 1. Flow of the Dirac eigenvalues for a fixed gauge field as a function of fj, for /3 = 1. 
(i) Two imaginary eigenvalues (and their partners) merge and move into the complex plane to 
form a quadruplet, (ii) A pair of imaginary eigenvalues merges at zero and becomes a pair of real 
eigenvalues, (iii) Two real eigenvalues (and their partners) merge and move into the complex plane 
to form a quadruplet. 



For j3 = 4 it follows from the symmetry (2.16b) that the nonzero eigenvalues come 
in quadruplets A, —A, A*,— A* (with eigenstates ip,j5ip,C^5tp*,C^p*) if A is complex and 
that the Kramers degeneracy is removed once /i is switched on. Hence we do not need the 
mechanism of figure l(i) and the eigenvalues move into the complex plane immediately. As 
in the case of /3 = 2, exactly real or imaginary eigenvalues only occur accidentally. The 
topological zero modes are twofold degenerate for both [x = and /x / 0. They change 
smoothly with \i and are eigenstates of 75 with the same chirality. 

3.2 Singular values 

Let us now consider the operator D^D. Since this operator is Hermitian with nonnegative 
eigenvalues, we write its eigenvalue equation in the form 

D^Dip n = tl^n ■ (3.4) 

The£ n are real and nonnegative. They are called the singular values of D, the name coming 
from the singular value decomposition of a non-Hermitian matrix. 20 The operators D^D 
and DD^ share all nonzero eigenvalues, since (3.4) implies 

Drf(D<p n ) = £{D<p n ), (3.5) 

and similarly the other way around. 

At \i = the A n and £ ra are trivially related by £ n = |A n |, and therefore the nonzero 
singular values are twofold degenerate (for (3 = 1 and 2) or fourfold degenerate (for f3 = 4). 
At fi 7^ there is no simple relation between the A n and £ n , i.e., knowing only the A n we 
cannot compute the £ n , and vice versa. As soon as /i is turned on the twofold degeneracy 
(for (3 = 1 and 2) is removed, and the fourfold degeneracy (for f3 = 4) is reduced to a 
twofold degeneracy, see the argument after (B.33). 

20 In this paper we assume that the extension of the singular value decomposition to non-Hermitian 
operators is straightforward and skip the mathematical foundations. 
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Figure 2. Schematic flow of the singular values of D for a fixed gauge field as a function of /i for 
the three symmetry classes. The numbers stand for the degeneracy of each singular value. 



The operator D^D commutes with 75, and therefore the states ip n have definite chiral- 
ity, 75</? n = i^n (or can be so arranged if the singular values are degenerate). Now assume 
£ n > and define (p n = ^~ 1 D(/9 n , for which 

DD^<f n = ^DtfDtpn = £ n D<p n = tilvn , (3.6) 

l<o¥n = CSs^Vn = -^~ X Dj S (p n = =flf n , (3.7) 

i.e., (p n is an eigenstate of DD^ with the same eigenvalue but chirality opposite to that 
of ip n . Therefore the number of right-handed (left-handed) nonzero modes of D^D equals 
the number of left-handed (right-handed) nonzero modes of DD^ . Since these operators 
coincide at fi = 0, the number of right-handed and left-handed nonzero modes of D^D 
is equal at \x = 0. As fj, is turned on from zero, a right-handed (or left-handed) nonzero 
mode of D^D changes its form smoothly but stays right-handed (or left-handed) since the 
eigenvalue of 75 is discrete. Therefore the numbers of right-handed and left-handed nonzero 
modes of D^D and DD^ are all equal (assuming that there are no accidental zero modes). 

If D has an eigenvalue equal to zero, Dip = trivially implies D^Dip = 0, i.e., there is 
a corresponding singular value of D equal to zero, and the zero mode of D*D is the same 
as that of D. If ip is not a zero mode of D, it cannot be a zero mode of D'D either since 
(ip\D< D\ip) 7^ 0. As discussed in section 3.1, all zero modes of D and therefore of D^D 
generically have the same chirality. Using D{^ = —D{—ji) and our earlier observation 
that the number of zero modes is stable as a function of fi, we conclude that (i) the operator 
DDt has the same number of zero modes as D^D, (ii) the zero modes of DD^ are equal to 
those of D(—fi), and (hi) the chirality of the zero modes of D^D and DD^ is the same. 21 

In figure 2 we schematically illustrate the singular value flow in QCD-like theories with 
(5 = 1, 2, and 4. All eigenstates of D^D have definite chirality, as indicated by R and L in 
the figure. For /3 = 4 the zero modes of D, and therefore of D^D, are twofold degenerate 
for both // = and \x / 0. 

21 Conclusion (iii) is not necessarily valid if there are accidental zero modes, while (i) and (ii) are always 
valid, see appendix D. 
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3.3 Dirac operator with chiral chemical potential 

In this subsection we derive a relation between the singular values and the eigenvalues of 
the Dirac operator with chiral (or axial) chemical potential defined by 

#5 Of) = luD v + /U7475 • (3.8) 

This operator was introduced in the context of the chiral magnetic effect [58]. 

Since D^{jj) is anti-Hermitian its eigenvalues are purely imaginary, and since it anti- 
commutes with 75 its nonzero eigenvalues come in pairs with opposite sign. Now, for any 
right-handed (left-handed) spinor ip R (<Pl) we have 

Ds(fj<) 2 <PR = {ptvDv + m<ah){ivD v + ^7475)^ 

= {luD u + /i7475)(7^A, + ^i)^R 

= -D(ji?D(jj,)<p R (3.9) 

and similarly 

D 5 (fj,) 2 ip L = -D(ji)D(ji?<p L , (3.10) 

D s (-n) 2 <p R = -D(jj,)D(ny<p R , (3.11) 

D B (-fj) 2 (p L = -D(ji?D{n)tp L , (3.12) 

which can also be expressed in terms of the decompositions 

D b (nf = -D^DPr - DD^P L , (3.13) 

M-V? = -D^DP L - DD^P R , (3.14) 

or equivalently 

D^D = -D 5 (fi) 2 P R -D 5 (-v) 2 P L: (3.15) 

DD* = -D 5 {n) 2 P L - D 5 (-fj,) 2 P R . (3.16) 

Let us denote the singular values corresponding to the right-handed (left-handed) nonzero 
modes of D^D by £, Rn (£,Ln)- From the above arguments we conclude that the nonzero 
eigenvalues of D^(fi) and D^{—^i) are given by the sets {±i£ Rn } and {±i£i, n }, respectively. 

Disregarding accidental zero modes, the zero modes of D^D and DD< have only one 
chirality. If they are all right-handed, the zero modes of D^{fj) (D5 (—//)) are right-handed 
and given by those of D^D (DD*). If they are all left-handed, the zero modes of Ds(fj,) 
(Ds(— //)) are left-handed and given by those of DD^ (D^D). 

Since Z?s(/x)^ = —D B (/j) and {D^(fj 1 ),^} = 0, the Banks-Casher relation for the chiral 
condensate can be extended to nonzero chiral chemical potential without difficulty, 

<V^> = 717*5(0) , (3.17) 

where /35(A) is the spectral density of D^(fi) with fermionic weight det JV f D^{fj). 
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4 Banks-Casher-type relations 



In this section we derive Banks-Casher-type relations for the three different theories with 
(3 = 1, 2, and 4. In each particular condensate is related to the density of the 

singular values at the origin. For our derivations to be correct it is important that in the 
computation of the partition function (including source term for the desired condensate) 
the fermionic measure is positive definite. If this requirement is not met a probabilistic 
interpretation is not possible and a number of subtleties arise that are discussed in some 
detail in appendix C. In the present section we restrict ourselves to cases in which the 
measure is positive definite. 

4.1 Two-color QCD {(3 = 1) 

For simplicity we take Jr = jrI and Jl = JlI with / given in (A. 9) and numbers jj?, ji 
that for the time being can be complex. It follows from (2.10) and the discussion after 
(2.12) that the measure is positive definite only if the combination e %e {—jRj* L ) N ^ 2 is real 
and positive. If we assume real sources and = 0, this condition is always satisfied for 
jn = —jl (corresponding to the scalar diquark condensate). If Nf/2 is odd, it is violated 
f° r 3R = 3l (corresponding to the pseudoscalar diquark condensate). In the following we 
therefore set jr = —ji = j with j real, 22 which implies J^JlPl + JrJrPr = j 2 ^-N f - From 
(2.10) we then obtain 

m = (de,<vw»+i 2 )) YM = (IF^T" 2 )™ 

n 

Introducing the notation 

(0) ^ZU)( 0de, "" /2(D,D + i2) )vM < 4 ' 2 ' 

for expectation values in the presence of a diquark source, we define the density of the 
nonzero 23 singular values of the Dirac operator by 

p sv (0= hm ^(y; *(£-&»)) for £>0, (4.3) 

Vi-s>oo I/4 \ I j=0 

where V4 is the space-time volume. The scalar diquark condensate can then be expressed 
in terms of this density at the origin, 

1 d 

(t/j T C-f 5 T 2 Iip) = lim lim — — lnZ(j) 



1 N f 2j 
lim lim ( > -s — 



Nf f'°° 2i 
J 1 dCp sv (0 .lim. 



3 



/-+o+ i 2 +j 2 



^vrp sv (0) , (4.4) 



22 We will comment on the case ju — jl below. 

23 See appendix C.2 for a discussion of the zero-mode contribution. 
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which is similar to the Banks-Casher relation for the chiral condensate in terms of the 
Dirac eigenvalue density at zero. Our analysis can be extended to nonzero quark masses. 
Assuming for simplicity M = ml with real to, we again obtain (4.1), but with D and D* 
replaced by D + m and D> + m. This means that (4.4) continues to hold, except that the 
singular values are now those of D + m. 

Some comments are in order. First, (4.4) holds at fx = and /j / 0, whereas the 
original Banks-Casher relation only holds at fi = 0. Second, in the derivation of (4.4) we 
have tacitly dropped the contribution of the singular values equal to zero. As discussed in 
appendix C, this is only justified if the measure is positive definite. Third, the integral in 
the third line of (4.4) needs a proper UV regularization. This was discussed carefully for the 
original Banks-Casher relation in [4] and works in exactly the same way here. The point is 
that the UV-divergent part disappears in the limit j — > + . Fourth, we observe that setting 
3R = 3l = j we would have obtained det W = det(D^ D+j 2 ) just as in (4.1), which formally 
would have led to (4.4) but with the pseudoscalar diquark condensate CT2I1P) on the 
l.h.s. instead. However, for odd Nf/2 the measure is not positive definite for jr = ji, which 
invalidates the Banks-Casher relation for the pseudoscalar condensate, see appendix C.2. 
For even Nf/2 the measure is positive for both types of sources, and what condensate is 
given by p sv (0) depends on the choice of sources we add. 24 In appendix C.3 we show that 
this is not in contradiction with QCD inequalities. Finally, we note that for two flavors 
and in the chiral limit, the relation (4.4) was obtained earlier by Fukushima [49], see also 
[59, 60]. Our result differs from that of [49] by a factor of 1/2. The contribution of the 
zero modes and the positive definiteness of the measure were not addressed in [49] . 

At fj, = 0, the spectra of the eigenvalues and singular values are identical, and so are 
the spectral densities at the origin. This implies that the chiral condensate (for m — > 0) 
and the diquark condensate (for j — > 0) are of the same magnitude, which is consistent 
with the fact that under a global SU(2Nf) rotation these condensates can be rotated into 
each other. As is well known, /1 ^ breaks this degeneracy. As /x increases (for to 7^ 0), 
the diquark condensate remains exactly zero until a critical value /i c = m n /2 is reached, 
and then starts growing for fi > \i c [20]. This behavior can be naturally understood by 
our Banks-Casher-type relation. For m 7^ the relation reads (iftip) oc p sv (0;m), where 
Psv(A;to) stands for the singular value density of D{p) + m, as remarked below (4.4). For 
sufficiently small fx, all eigenvalues of -D(/i) are still localized near the imaginary axis, and 
the density of the near-zero modes of D(fi) + to is zero. As p increases, the eigenvalues 
spread out more, and for ll > \i c there is a nonzero density of eigenvalues at ±m. This 
signals a nonzero p sv (0;m), i.e., the onset of diquark condensation. 25 

It is also possible to express the partition function in terms of the Dirac operator 
with chiral chemical potential D§{p) defined in (3.8). This is most easily shown working 

24 More precisely, the magnitude of the condensate is given by p sv (0) and its orientation by the external 
sources. A simple analog is the Ising model, where the direction of the spontaneous magnetization at zero 
temperature depends on the direction of the (infinitesimal) external magnetic field. See also appendix C. 

25 A similar discussion starting from a different method can be found in [61]. 
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backwards starting from (4.1) 



m = (u& +f) Nf/2 ) YM = {j nm [IT(&, +i 2 )II / (^+i 



2\ 



N f /2 



YM 



det^^s^+jjdet^^sC-^+i))^, (4-5) 

where is the number of topological zero modes of D^D, the primed products are only 
over nonzero singular values, and in the last line we have used the relationships between 
the eigenvalues of D$ and the singular values derived in section 3.3. 

4.2 QCD with isospin chemical potential {(3 = 2) 

It follows from (2.14) that the measure is positive definite only if the combination e l8 (—p\*) 
is real and positive. We therefore choose 6 = and p = —A = j with real j (another choice 
for which the measure is not positive definite is considered in appendix C.l). The partition 
function is then 

Z(j) = (det(D*D + j 2 )) YU , (4.6) 

and by a calculation analogous to section 4.1 we obtain the pion condensate 

(uj5d - dj 5 u) = 7rp sv (0) . (4.7) 

Similar comments as in section 4.1 apply. In particular, we have dropped the contributions 
of the zero modes (which is justified because the measure is positive definite), and a proper 
UV regularization is understood. As in section 4.1 we can also express the partition function 
in terms of D^(p) and obtain the same result as in (4.5). 

4.3 QCD with adjoint fermions (/3 = 4) 

We now take Jr = j#l and Jl = j'lI. The measure in (2.20) is then positive definite 
only if the combination e t9 ^ Nc (—jnj^ j ) N f is real and positive. As in section 4.1 we therefore 
assume 6 = and set jr = —ji = j with real j (see appendix C.l for another choice). The 
partition function is then 

Z( J ) = (det^/ 2 ( J Dt jD + i 2 ) \ ( 4. 8) 

\ / YM 

and a calculation analogous to section 4.1 leads to 

(^ T C 75 ^> = -firp BV (0) . (4.9) 

Again, similar comments as in section 4.1 apply, and we could also have expressed the 
partition function as in (4.5). 

5 Chiral Lagrangians with diquark source (f3 = 1) 

From now on we concentrate on two-color QCD (/3 = 1) for simplicity and brevity. Results 
similar to those obtained in sections 5-7 could also be derived for the theories with (3 = 2 
and f3 = 4 using the methods employed here, but we will not pursue this here. 
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5.1 Three different regimes: low, intermediate, and high density 

We now construct the low-energy effective theory for two-color QCD at nonzero density in 
the presence of a diquark source and in the chiral limit. There are actually three different 
effective theories, applicable at low, intermediate, and high density. In the following we 
will refer to them as L, I, and H, respectively. The effective theory L is constructed under 
the assumption of maximal chiral symmetry breaking at low density [19, 20], while the 
effective theory I is constructed assuming the conjectured BEC-BCS crossover discussed 
in the introduction. 26 On the other hand, the effective theory H is constructed based on a 
rigorous weak-coupling analysis at high density [15]. 

In this subsection we discuss how the three density regimes differ in their patterns of 
spontaneous symmetry breaking and the number of Nambu-Goldstone (NG) modes and 
comment on the connection between the three regimes. The effective theories themselves 
and mass formulas for the NG modes will then be derived in the next three subsections, 
and their domains of validity will be discussed in section 5.5. 

At very low density one can start from the pattern of chiral symmetry breaking at zero 
density, 

SU(2N f ) Sp(2N f ) , (5.1) 

and treat the chemical potential and the diquark source as a small perturbation. This is 
the approach taken in [19, 20]. The NG modes of the theory are the same as those of 
the zero-density theory, i.e., they are collected in a field S that parametrizes the coset 
space SU(2JV»/Sp(2JV». SU(iV) and Sp(iV) have N 2 - 1 and N(N + l)/2 generators, 
respectively, hence the number of NG modes in this regime is Nj(2Nf — 1) — 1. One should 
keep in mind, however, that some of these modes acquire a mass as \i is increased. 

At intermediate density, when \i can no longer be treated as a small perturbation, 
we first recall that \x breaks the original SU(2iVf) symmetry to SU (NAl x SU (Nf)n x 
\J(1)b- A diquark condensate then breaks this symmetry to Sp(Nf)i, x Sp(Nf)n so that 
the symmetry-breaking pattern is now 

SU(JV» L x SU(JV» fl x U(l) fl ->• Sp(JV»£ x Sp(N f ) R . (5.2) 

The corresponding NG modes are S£,S^ 6 S\J(Nf)/Sp(Nj) and V G U(l), and the total 
number of NG modes in this regime is Nf(Nf — 1) — 1. 

At very high density the U(1) J 4 anomaly is suppressed due to the screening of instan- 
tons [66-68]. We therefore need to take the original U(1)a symmetry of the action into 
account. It is no longer broken explicitly by the anomaly but spontaneously by the diquark 
condensate so that the symmetry-breaking pattern is 

SU(iV» L x SV(N f ) R x U(l)s x XJ(1) A Sp(N f ) L x Sp(N f ) R . (5.3) 

26 From this point of view, our results from the effective theory I below should be viewed as predictions to 
be verified in future lattice simulations, which would then (dis)confirm the conjectured BEC-BCS crossover. 
Other possibilities for effective theories in the intermediate-density regime have been considered earlier, see, 
e.g., [62, 63]. In the model of [62] certain vector mesons could become massless at nonzero density, while 
in [63] this does not happen. Lattice studies on this issue are inconclusive [64, 65]. Here, we assume that 
all vector mesons remain massive at all densities. 
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Figure 3. Schematic /^-dependence of the masses of the NG modes. See text for details. 



The NG modes are the same as at intermediate density, except that there is an additional 
NG mode corresponding to U(1)a which we call rj' . In other words, the rj' mass has become 
negligible. Hence the total number of NG modes in this regime is Nf(Nf — 1). 

Let us make some qualitative comments on the connection between the three regimes, 
starting at zero density. Assuming the diquark sources to be infmitesimally small, there 
are Nf(2Nf — 1) — 1 massless NG modes at fi = 0. As n is increased, Nf(Nf — 1) — 1 
of them remain massless while Nj of them acquire a //-dependent mass [20]. Starting 
from the effective theory at zero density, these Nj modes can be integrated out. This 
yields the effective theory at intermediate density, which does not depend explicitly on /z 
but contains parameters (i.e., low-energy constants) that have acquired a //-dependence 
through the integrating-out of the massive modes. A similar argument applies starting at 
high density. As fi is lowered, the U(l)yi anomaly reappears so that the rf becomes massive 
and can be integrated out. These comments are illustrated in figure 3 and will be made 
more quantitative in the next subsections. 

Before proceeding, let us discuss a somewhat subtle issue regarding the parity of the 
diquark condensate. In general, when constructing an effective theory, one starts from 
an assumption of how the symmetries of the theory are broken, i.e., this assumption is 
an input to the effective theory. In the construction of the effective theories below, one 
assumption we have to put in is whether the scalar or the pseudoscalar diquark condensate 
minimizes the ground-state energy. For nonzero mass and zero diquark source it was 
shown by QCD inequalities that if a diquark condensate forms, it does so in the scalar 
channel [19]. However, for zero mass and nonzero diquark source QCD inequalities do 
not provide any information on this question, see appendix C.3. The assumption we will 
make is that the diquark condensate again forms in the scalar channel. This assumption 
is based on instanton dynamics at high density, see, e.g., the review [69]. The instanton 
vertex is cdet(^V-R) + h- c - with c > 0, where we have Nf legs each for i(jr and 
[70]. Taking the expectation value with respect to the diquark-condensed ground state, 
the contribution of the instanton vertex to the energy is d (tp^ L ^ L ) Nf ^ 2 (i/}riPr) N ^ 2 with 
d > 0. For Nf = An + 2 with n E N we obtain c'((i)\^ ] L ) (^r^r)) 2 " 1 x (V'lV'i) (tpRipR), 
which is negative for (iPl^Pl) = — {ipRipR.) but positive for (i/jliPl) = {ipRipR.)- Therefore, 
the positive-parity state is favored by instantons. We assume that this argument carries 
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over to low density based on the conjectured BEC-BCS crossover. For Nf = An the 
argument does not apply since the contribution to the energy is cf((x/j L tp L ) (ipRipR)) 2 '' 1 , 
which is not affected by a relative sign between (V'lV'I/) an d {^r^r)- Although in the 
construction of the effective theories below we do not distinguish between Nf = An and 
Nf = An + 2, we will see that for Nf = An the positive- and negative-parity states turn 
out to be degenerate, while for Nf = An + 2 the pseudoscalar condensate is suppressed or 
vanishes completely. This is also consistent with our analysis in section 4.1: For Nf = An 
the microscopic theory always has a positive definite measure so that we can derive Banks- 
Casher-type relations for both the scalar and the pseudoscalar diquark condensate, showing 
that their magnitudes are equal. For Nf = An + 2 the pseudoscalar diquark source leads to 
an indefinite measure in the microscopic theory, and therefore we cannot conclude anything 
about the relative magnitude of the two condensates. 

5.2 Effective theory L at low density 

At low density the effective theory derived in [19, 20] uses as degrees of freedom the NG 
modes corresponding to the symmetry-breaking pattern (5.1) at zero density. We briefly 
review and extend the relevant results here. In the chiral limit, the leading-order effective 
Lagrangian is given by 

F 2 

£eff = -y tr(V„£ V„Et) _ $ L Re tr ( JE) (5.4) 



with 



V v ^ = d v E-nS vQ {BE + EB), (5.5) 
V V T) = 8 U ^ + fi5 u0 (^B + B^), (5.6) 



The field E parametrizes the coset space SU(2Nf) /Sp(2iVy), 



T , , „ , j u \ i'iTT a T° 




S C^iU' with E,,= ( m , ) and U = exp [-fir) » (5.8) 

where the T a (a = 1, . . . , N f (2N f - 1) - 1) are the generators of SU(27V / )/Sp(2A r / ). They 
are Hermitian and satisfy tr(T a ) = and tv(T a T b ) = 5^. A sum over the repeated index 
a in (5.8) is understood. 

There are four minor differences with respect to [20]. First, only the case of Jr = —Jl 
was considered there. Second, our <3?l corresponds to their G. Third, the analysis in 
[20] also includes a nonzero quark mass m. We only consider m = 0, in which case the 
chiral condensate immediately disappears as /i is switched on, while the diquark condensate 
immediately assumes its full value, see section 12 of [20]. Fourth, to be consistent with 
the rest of the current paper, our convention for the diquark source and consequently 
for Erf differs from [20]. This follows from requiring parity invariance of the microscopic 
Lagrangian (which implies Jl — Jr under parity), of the diquark condensate (which 



- 20 - 



determines the form of S^), and of the effective theory (which determines the form of J 
and the second term in (5.4)). As a consistency check, we note that the minimum of the 
action is obtained for £ = and Jr = —Jl = jl with j real and positive. 

In (5.4) there are two low-energy constants. F is the common decay constant of all 
NG modes, and the positive parameter <I>l is the magnitude of the diquark condensate per 
flavor and handedness in the absence of sources and at fi = 0, 

<*>L = ^ MCt 2 I^)\-j =0 ^ =0 (i = L,R). (5.9) 

It is important to note that F and <I>l do not depend on ft. 

For Jr = — Jl = jl the masses of the NG modes have been computed in [20, eq. (101)], 
and in the chiral limit we have (p = a = ir/2 and = Vj3>/ ' F in these expressions. Hence 
there are two types of NG modes, 

type 1: mass = ^Jj^jF 2 (Nf - N f - 1 modes) , (5.10a) 

type 2: mass = \Jj^lF 2 + (2fi) 2 (Nf modes) . (5.10b) 

While the type-1 modes are massless in the j — > limit, (5.10b) shows that ft is an explicit 
symmetry-breaking parameter that makes the type-2 modes massive even for j = 0. 

Let us comment on the source Jl = Jr for the pseudoscalar diquark condensate. 27 For 
Nf = An we define u = diag(il4 n , l4 n ). Since u £ SU(8n) we can redefine U — > ull in (5.8). 
The measure is not changed by this transformation, but in the term JE in (5.4) the sign 
of Jl in (5.7) is flipped. This means that the partition functions (and hence the energies) 
for Jl = — Jr and for Jl = Jr are exactly equal. This is consistent with the observation 
that the instanton vertex does not prefer one condensate over the other for Nf = An, see 
the discussion at the end of section 5.1. For Nf = 2 we find that a pseudoscalar diquark 
source term drops out from (5.4) since Retr(JE) = Retr[diag(J, I)C/diag(I, — I)U T ] = 
for U € SU(4). 28 Hence the pseudoscalar diquark condensate is zero (in this order of 
the low-energy expansion), which is again consistent with the instanton-based argument 
in section 5.1. For Nf = An + 2 with n > 1 the situation is more complicated. While 
we currently cannot make any definite analytical statement, numerical experimentation 
indicates that for Jr = Jl = jl the minimum of the energy is larger than for Jr = 
—Jl= jl- If true, this would mean that by changing the diquark source we can generate 
a nonzero pseudoscalar diquark condensate whose magnitude is smaller than that of the 
scalar condensate. 

5.3 Effective theory H at high density 

For technical reasons we now proceed to the regime of very high density. The effective chiral 
Lagrangian including mass term for this case was derived in [15]. Here, we are interested 
in diquark sources and therefore set the mass term to zero. 

27 This is equivalent to Jl = —Jr with 8 = Nj-ty/2, see section 2.1, i.e., instead of studying the J- 
dependence we could equivalently study the ^-dependence. This statement applies at any density. 

28 This can be shown using an explicit parameterization of the coset space SU(4)/Sp(4), see, e.g., [71]. 
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We first consider the case of Nf > 4. The symmetry-breaking pattern is given in (5.3). 
By forming linear combinations of the generators of V(1)b and U(1)a we can switch from 
U(1)b x U(1)a to U(1)l x XJ(1)r. We parametrize the NG modes as 



Si = UilUf with Ui = exp f ) (i = L,R), (5.11 




/- exp [ v a'/J • " ,,x,, l • (5 - l2) 



The parameterization of the Si is similar to (5.8), except that the T a (a = 1, . . . , Nf(Nf — 
l)/2 — 1) are now the Hermitian generators of SU(./Vy)/Sp(iVj), again satisfying tr(T a ) = 
and tr(T a T b ) = 5 ab 29 Using U l l\jJ = Ufl [20] we can also write Si = Ufl. 
The quarks transform under S\J(N f ) L x SXJ(N f ) R x U(1) L x U(1)r as 

i>L^e iaL g L ^ L , ^ R ^j aR g R ^R, (5.13) 

where gi 6 SU(iVy-)i and e iQ!i £ U(l)j (i = L,R). The NG modes therefore transform as 

Si -> 5iS i5 f (i = L,R), L — >• Le 2iQi , i? -»• i?e 2iQfl . (5.14) 

The transformation properties of and J/j are determined by requiring that the Lagran- 
gian (2.3) be invariant under the flavor symmetries. This implies 

Jl 9lJ L g{ e~ 2iaL , Jr -> <&Jjtf J, e" 2ia * . (5.15) 

Therefore the invariant real combination linear in Jl and Jr is uniquely determined to be 

Re [Ltr(J L S L ) - R^JrEr)] . (5.16) 

As for theory L, we required parity invariance of the microscopic theory (implying Jl 
—Jr), of the diquark condensate (implying S^ S# and L o i£), and of the effective the- 
ory (leading to the relative factor of —1 in (5.16)). The leading-order effective Lagrangian 
in the presence of diquark sources and in the chiral limit is thus given by 



*-eff 



^(|9 L| 2 + vl\dM 2 ) + y tr (|3 S L | 2 + ^S L | 2 ) + (L O fl) 
- $ H Retr(J L LS L - J flj RS R ) - ^m 2 lst Re (L^// 2 , (5.17) 

where /o, / and vo, v are low-energy constants that correspond to the decay constants and 
velocities of the NG modes, respectively. The latter are generally different from unity (i.e., 
the speed of light) since Lorentz invariance is lost at /i 7^ 0. The minus sign in front of the 
positive low-energy constant $h is chosen so that the minimum of the action is obtained 
for L = R = 1, Til = Er = I, and Jr = —Jl = jl with j real and positive. Note that all 
29 The T a are related to the X a in [15] by T a = X a /JWf. 
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low-energy "constants" in depend on /i. Their relation to physical observables will be 
discussed in section 5.6. 

In (5.17) we have also included a term that corresponds to the single-instanton contri- 
bution to the r/ mass, parametrized by mi ns t. This term is symmetric under the anomaly- 
free subgroup (Z^a^ )a of U(1)a- From the symmetry point of view all terms of the form 
Re(tfR) nN f/ 2 (n > 1) are allowed and contribute to the 7/ mass. Microscopically, these 
terms correspond to n-instanton vertices. At sufficiently large the instanton ensemble 
can be regarded as a dilute gas [66-68] which does not form instanton molecules [72]. The 
diluteness of the instanton gas is parametrized by the dimensionless quantity (proportional 
to the instanton density) a oc (A QC D/n) b{Nf) < 1 with b(N f ) = (22 - 2N f )/3. Since the 
probability ~ a 11 that n instantons (n > 2) are at the same point is highly suppressed 
we can neglect the multi-instanton vertices and only keep the one-instanton contribution. 
Note that mi ns t is a decreasing function of /i, with mi nst — > for fj, — > 00 [66-68]. 

It is convenient to combine the U(l)j field with S\J(Nf)i/Sp(Nf)i by defining 

/ i7r A T A\ 

S L = LS L = exp — k / and (L O R) , (5.18) 



fA 

where the T A (A = 0, . . . , N f (N f -l)/2-l) are now the generators of U(iV/)/Sp(iV/) with 
T° = t/y/Nf so that ti(T A T B ) = 5 A b 30 We also defined f A = f for A > 1. To second 
order in the 7r-fields we have 

/ Tr A 7r B T A T B \ 
ReE i= 1-^yL I. (5.19) 

Assuming Jr = —Jl = 3 1 with j real and positive, this yields a Gell-Mann-Oakes-Renner 
(GOR) type mass formula for the tt a . As in theory L there are two types of NG modes, 



typel: m A = Jj$ H /f% (Nj - N f - 1 modes) , (5.20a) 



type 2: rry = yjj$n/f$ + mf nst (1 mode) . (5.20b) 

Note that there are two type-1 modes for each A > 1, but only a single one for A = 0. 
Note also the similarity with (5.10). Now m; nst plays the role of the symmetry-breaking 
parameter which makes the type-2 mode massive as \x is lowered. 

Let us now consider the case of Nf = 2, in which and X# are absent because 
SU(2) ~ Sp(2). Thus the effective Lagrangian contains only the fields L and R. Since 
any 2x2 antisymmetric matrix is proportional to I we can write Jl = jil and Jr = 
JrI (Jl, 3r £ C) without loss of generality. These terms transform as 

JL -> 3l e~ 2iaL , 3R -> OR e~ 2iaR , (5.21) 

which follows from (5.15) and g*Igj = (det g*)I = I (i = L,R). Therefore the invariant 
real combination linear in jr and jl is given by 

Re(j L L-j R R). (5.22) 



30 We use the convention that uppercase indices corresponding to U(iV/ )/Sp(JV/) start at zero, while 
lowercase indices corresponding to SU(iV/)/Sp(iVf ) start at 1. 
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Using the same parameterization of L and R as in (5.12), the effective Lagrangian for 
Nf = 2 reads 

£rfE = /o [\doL\ 2 + vl\8iL\ 2 + (L o fl)] + 2d> H Re(j L ^ - Jflfl) - ftmLt Re(L^R) , 

(5.23) 

where the plus sign in front of and the factor of —1 between the two terms following 
it have been chosen so that the minimum of the action is obtained for L = R = 1 and 
jr = —ji = j with j real and positive. The GOR-type relation for this case is identical to 
(5.20) with A = 0. 

We again comment on the case of Jl = Jr- Note first that for mi nst = (i.e., at 
asymptotically high density) the fields L and R in (5.17) or (5.23) can rotate indepen- 
dently, and hence the left and right diquark condensates in the ground state can be rotated 
separately by varying the directions of Jl and Jr. This is no longer true when mi ns t ^ 0. 
Since the anomaly term favors L = R energetically, the fields L and R can no longer rotate 
independently. Now let us again discuss the various cases of Nf. For Nf = 4n we can 
redefine Ul — > uUl in (5.11) with u = diag(i]U n ) G SU(4n), which flips the sign of El and 
thus absorbs a sign flip of Jl in (5.17). So again the energies for Jl = —Jr and Jl = Jr 
are equal, in agreement with the instanton-based argument in section 5.1. For Nf = 2 
and jr = —jl > the diquark source term and the anomaly term in (5.23) are minimized 
simultaneously at L = R = 1, and thus the ground state is not changed by the anomaly 
term. However, for Jr = j'l > there is a competition between these two terms, and they 
cannot be minimized simultaneously. Therefore the pseudoscalar diquark condensate can 
be realized only if the diquark sources are strong enough to overcome the penalty due to 
the anomaly term. For Nf = 4n + 2 with n > 1 the situation is similar but a bit more 
complicated. For Jr = —Jl = jl both the diquark source term and the anomaly term 
in (5.17) can be minimized simultaneously. For Jr = Jl = jl we suspect, although we 
currently cannot prove it analytically, that this cannot be done. A more quantitative study 
is required to determine the magnitude of the pseudoscalar diquark condensate in this case. 

5.4 Effective theory I at intermediate density 

At intermediate density, the coupling constant is not small enough to treat instantons as 
a dilute screened gas, and hence the \J(1)a anomaly can no longer be treated as a small 
perturbation. In other words, mi nst increases as // is lowered so that the r/ mass in (5.20b) 
is not necessarily small, implying that the r( should be integrated out from the effective 
theory. Technically, this means that L and R should be replaced by a single U(l) field V, 
i.e., the effective Lagrangian (5.17) for Nf > 4 changes to 

4fl = Wo [\d V\ 2 + v 2 \diV\ 2 ] + £ tr [\d Z L \ 2 + v 2 \diZ L \ 2 + (L O R)} 

- $i Re [V tr( J L S L - JrZr)} (5.24) 

corresponding to the symmetry-breaking pattern (5.2). We have renamed the low-energy 
"constants" to distinguish them from those of theory H. Note that they again depend on 
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\i. In section 5.6 we will discuss how they are related to the low-energy constants of theory 
L and H at low and high density, respectively. 

The NG modes ^l/r are parametrized as in (5.11) with / replaced by /, while 

V = exp | jjL 1 . (5.25) 
To second order in the 7r-fields we now have 



7T 



\2 ^0 ^arpa 



vrP / 7rfT a ir?i$T a T b 



Re^iV) = I 1 - - -LKJ^ " I "^; 9 - I /. (5.26) 

Assuming again J/j = — Jf, = jl with j real and positive, we obtain a GOR-type mass 
formula analogous to (5.20a), 



m A = \Jj$l/f A (Nf -N f -1 modes) , (5.27) 

where Ja = / for A > 1. Note that we have only type-1 modes in theory I. 
For Nf = 2 the effective Lagrangian changes to 

£* ff = 2/ 2 [|<9 y| 2 + v 2 \diV\ 2 ] + 2*i Re[(j L - j^F] , (5.28) 

and the GOR-type relation for the single NG mode is tuq = f$ as in (5.27). 

Finally, we again consider the case of Jl = Jr. For Nf = 4n the argument and the 
conclusion are exactly the same as in theory H. For Nf = 2 the pseudoscalar diquark 
source term in (5.28) drops out trivially since Jl = jn, and hence the pseudoscalar diquark 
condensate is zero (in this order of the low-energy expansion) as in theory L. For Nf = 4n+2 
we suspect, although we currently cannot prove it analytically, that for Jr = Jl = jl the 
minimum of the ground-state energy is larger than for Jr = — Jl = jl- If true, we are led 
to the same conclusion as in theory L. 

5.5 Domains of validity 

In the following discussion we assume that the diquark sources are infinitesimal. In general, 
there are two conditions for an effective theory formulated in terms of NG modes to be 
applicable: (i) the masses of all NG modes must be much smaller than the mass scale of 
the lightest non-NG particle, and (ii) the typical scale p of observables computed within the 
effective theory must also be much smaller than m^. Of course, ni£ itself must be nonzero. 
To figure out the domains of validity of the three effective theories at nonzero density we 
must determine the mass scale mi of each theory, which generically is a function of [i. 

For the effective theory L we have m^(L) ~ A, where A is the mass of the lightest 
non-NG particle at zero density. As [i is increased from zero, some of the NG modes of L 
acquire a mass proportional to fx. The effective theory I is obtained from L by integrating 
out these modes so that for I at low density we have me(I) ~ fx. 

The situation at high density is somewhat more complicated. At asymptotically high 
density the r/ is massless and A ^ ll. For the effective theory H we have m^(H) ~ A. 
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effective theory 


condition (i) 


condition (ii) 


Li 


/i<A 


p <C A 


I (scenario 1) 




p <^ /j, ior ^ a 

p < A for ^bcs < /" < /"c 

p < m,,/ for ;U > ^ c 


I (scenario 2) 




p <ZC h for /i < A 
p < my for /i > /i B cs 


H 




p < A 



Table 1. Domains of validity of the three effective theories L, I, and H, see text for details. 

There are two ways to see this. First, A plays the role of a constituent quark mass so that 
the lightest non-NG particles (color singlet diquarks and mesons) weigh about 2A [22]. 
Second, the higher-order vertices in the effective Lagrangian are suppressed by 1/A, while 
loop integrals are suppressed by 1//a [73]. Since /a ~ H [74] and A <C the cutoff is A. 
Let us now lower the density so that the r/ becomes massive, but let \x be large enough so 
that we still have my <C A <C [66-68]. The effective theory I is now obtained from H 
by integrating out the rf so that for I in this regime we have m^(I) ~ my. As the density 
is lowered further there are two possible scenarios: 

1. There could be a "critical" chemical potential \x c at which my = A and below which 
my > A. In that case m,£(I) ~ A for hbcs ^5 f 1 ^5 He where hbcs is the chemical 
potential above which we are in the BCS regime. 31 

2. We could have < A for all /j, > hbcs- Then m^(I) ~ m ri i for all \i > hbcs- 

Since we only know the functions A(/i) and m^^/i) at asymptotically large fi we do not 
know which of these two scenarios is correct. 32 This can only be decided by a full dynam- 
ical calculation, e.g., in lattice QCD. However, the two scenarios can be combined in the 
statement rn^(I) ~ min(A(^), m^'(/^)) for [i > hbcs- 33 

Our discussion is summarized in table 1 and figure 4. Three comments are in order. 
First, condition (i) is always satisfied for the effective theory I since the Nf(Nf — 1) — 1 
NG modes shown in figure 3 are always massless. In other words, I is applicable at any \i 
as long as the scale of the observable is sufficiently small. Second, as the "intermediate" 
density is increased from "low" to "high", the mass scale of the lightest non-NG mode 
changes from /j, to my . Finding the precise //-dependence of mi in the intermediate region 
again requires a full dynamical calculation. All we currently know are the two limits fi 

31 Although there is no phase transition between the BEC and BCS regimes, we can define hbcs as 
the chemical potential above which the minimum of the dispersion relation of the fermionic quasiparticles 
changes from p = to p 7^ 0, with p being the momentum [75]. 

32 In these asymptotic functions we always have m r y < A, but this does not tell us anything about the 
regime where fi is not asymptotically large. 

33 In the preceding arguments we have completely ignored gluons, even though they are lighter than the 
Vj at sufficiently large fi [68], since their interaction with NG modes is assumed to be negligibly small, see 
the discussion in [15]. 
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Figure 4. Domains of validity of the three effective theories L, I, and H in scenario 1 (left) and 
scenario 2 (right), see text for details. 

and m^i. Third, in the overlap regions of the different effective theories (i.e., the green and 
orange areas in figure 4) one has a choice of which theory to use, but this choice depends 
on the observable. (This is similar to the choice between SU(2) and SU(3) flavor chiral 
perturbation theory in QCD as a function of the strange quark mass [76, 77].) For example, 
in the regime /i <C A one could also use the effective theory I, but this would only work 
for observables with p <C [i, whereas the effective theory L could be used for observables 
with p <C A. Similarly, the effective theory I could also be used in the regime mw <C A. 
This would only work for observables with p <C m„/ , whereas the effective theory H could 
be used for observables with p <C A. While at first sight it may seem that one should 
always work with the effective theory that allows for a larger range of observables, it may 
be technically simpler to work with the effective theory I if one is only interested in an 
observable for which this effective theory is valid. 

5.6 Matching of the low-energy constants 

Let us now comment on the relation between low-energy constants and physical observables, 
and on the matching of the low-energy constants between the different effective theories. 
Note that there are also low-energy constants corresponding to higher-orders in the effective 
Lagrangians, which we have not shown explicitly. 

Let us start with theory I, because this is simplest as the diquark source is the only 
symmetry-breaking perturbation. In that case /o and / are equal to the physical decay 
constants at a given /j, in the limit J — > 0, and 

$I = W f M Ct ^)\j= (' = L,R) ■ (5.29) 

For nonzero J there will be corrections (similar to the chiral corrections in chiral per- 
turbation theory) due to which the low-energy constants will deviate from the physical 
quantities. Next we consider theory L. It has two symmetry- breaking perturbations, [i and 
J. The low-energy constants do not depend on these external parameters. Rather, they 
are equal to the decay constants and the diquark condensate at fi = J = 0, see (5.9). Now 
let us recall that the theories L and I have overlapping domains of validity at low density. 
All physical quantities should be independent of which effective theory we use. As long 
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as we work at any finite order of the low-energy expansion in theory L, the results thus 
obtained could be different from those of theory I. However, the discrepancy will disappear 
if we sum up the contributions of the heavier NG modes (with mass ~ /i) in theory L to 
all orders. 34 We expect, for any fixed fi <C A, the relation 

| (V'V') I J=o = $i = + (Corrections due to the propagation of type-2 modes) , (5.30) 

and likewise for F, /o, and /. The relations between the low-energy constants of theory L 
and I can be made more precise by explicitly integrating out the type-2 modes of theory 
L. In the course of this procedure, the (initially ^-independent) low-energy constants ac- 
quire a /i-dependence in much the same way as the low-energy constants of SU(2) chiral 
perturbation theory acquire a dependence on the strange quark mass m s when kaons are 
integrated out of SU(3) chiral perturbation theory [76, 77, 79, 80]. According to these cal- 
culations, the corrections to the low-energy constants at 0(p 2 ) (F and B in the standard 
notation) due to the integrating-out of kaons become arbitrarily small when m s gets small. 
If we assume that this finding persists in our present context, we expect 

lim f ,f = F and lim $ r = $ L . (5.31) 

However, we do not expect this smooth matching to extend to the low-energy constants 
of higher orders, because it is known that in SU(2) chiral perturbation theory they receive 
corrections of the form (l/m s ) n with n > and thus blow up as m s — > 0. Hence they cannot 
reduce to the low-energy constants of SU(3) chiral perturbation theory [76, 77, 79, 80]. 

The discussion for theory H and its matching with theory I proceeds analogously. In 
theory H, the role of ll in theory L is now played by rrii nst , which in turn is a function of 
\x. For sufficiently high density the domains of validity of theory H and I overlap, and for 
any fixed /i (provided that m.i nst <C A) we expect the relation 

\(tpip)\j = Q = $j = $ H + (Corrections due to the propagation of 77 ) (5.32) 

to hold (and likewise for /o, fo and /, /). Based on our discussion at low density we now 
expect the matching (at any fixed li) 

lim JoJ = fo, f and lim $ r = 3> H . (5.33) 

However, there is a subtlety specific to high density. For li — > 00 the low-energy "constants" 
are actually infinite since fo, f ~ ll [74] and $ ~ Li 2 A/g [81]. Therefore f , /, and $ H 
cannot be defined as constants at fj, = 00. Accordingly, theory H cannot be defined at 
[x = 00 in the same way as theory L could be defined at fj, = 0. 35 That is why we defined 

34 This is similar to what is encountered in weak-coupling perturbation theory. In this case observables 
(such as the cross section) depend on the renormalization scale at any finite order of perturbation theory, 
but this dependence decreases as we go to higher orders [78]. 

35 This is not meant as a negative statement since the asymptotic behavior of the low-energy constants in 
terms of /1 is known. Similar situations occur when considering the large- iV c limit or scattering processes 
in the limit of high energies. 
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theory H and its low-energy "constants" at fixed (and finite) fi. This implies that mi ns t 
is also fixed and cannot be considered as an independent symmetry-breaking parameter 
anymore. However, we can start from theory H at a given fi and formally integrate out 
the rf to obtain theory I at the same value of fj,. The low-energy "constants" of I thus 
acquire a dependence on m; ns t, and we can now formally send mi ns t — > 0, still at the same 
fixed /x. This is how (5.33) should be understood. 

At the end of this section, let us comment on possible extensions. We could have 
performed a more comprehensive analysis with nonzero diquark sources, which would give 
a small mass proportional to \fj to the NG modes. However, in this paper we are only 
interested in the limit j — > 0, and therefore we have not performed such an analysis. 36 
In principle one could also add explicit quark masses. This makes the analysis even more 
complicated since the coset space one should use to construct the effective theory now 
depends not only on /j, but also on the quark masses. From the arguments presented in 
this section it should be clear how to proceed, but we do not pursue this issue further. 

6 Smilga-Stern-type relations (f3 = 1) 

In [2] Smilga and Stern computed the slope of the density of Dirac eigenvalues at the origin 
in the QCD vacuum (/3 = 2) using effective-theory techniques. Their result was confirmed 
by partially quenched chiral perturbation theory for degenerate [82] and nondegenerate [83] 
masses. Later it was generalized to theories with /3 = 1 and 4 at /i = [84]. In this section 
we adapt the method of [2] to the singular values at \x ^ 0, i.e., we compute the slope of 
the singular value density of the Dirac operator in two-color QCD at nonzero \i (f3 = 1), 
using the effective theories constructed in section 5. We will obtain three different results 
at infinite, intermediate, and zero density, respectively. In section 6.4 we will discuss the 
relation between these results as a function of \x. Throughout this section we work in the 
chiral limit for simplicity. 

6.1 Infinite density 

For technical reasons we now start at infinite density so that we can set mi ns t = in (5.17). 
As in earlier sections we set Jr = —Jl = J, where J is an antisymmetric Nf x Nf matrix 
that has Nf(Nf — l)/2 independent components. We can decompose J as 

J = /^i^ A = i/ + /^j a t a , (6.1) 

A a 

where the t A are the generators of U (N j) j 'Sp(iV f) and the ja are real parameters with 
io = j \J Nf . Such a decomposition is possible since the dimension of XJ(Nf)/Sp(Nf) is 
Nf(Nf — l)/2 and thus equal to the number of degrees of freedom of J, and since It A is 
antisymmetric for all A (which follows from t A I = I(t A ) T [20]). As before, the sum over 
A starts at 0, while the sum over a starts at 1. The t A are identical to the T A defined 
below (5.18), but we denote them by a different symbol (in agreement with the notation 

36 But see section 6.4 in which we are forced to consider the case of small but nonzero j to understand an 
apparent discontinuity of a very particular observable. 
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of [84]) since they are used in a different context: The T A are used to parametrize the 
NG modes of the effective theory living in \J(N f ) L /Sp(N f ) L x U(N f ) R /Sp(N f ) R , while the 
t A are used to parametrize the source J, which exists already in the microscopic theory. 
The choice of U (Nf) /Sp(NA to parametrize J is natural since it is the space in which the 
diquark condensate aligns depending on the choice of J. 
For Nf > 4 we now consider the scalar susceptibility 



K ab (j)= lim L dja d 3b lnZ(J) 
V4— s-oo I/4 



(6.2) 

all j a =0 



which we will calculate both from the microscopic theory (two-color QCD) and from the 
low-energy effective theory. 

Let us start on the QCD side. Assuming that at high density there are no zero modes, 
(2.10) yields 

Z(J) = (det 1 / 2 ^ + JU)) ym = ( IIdet 1/2 (^ + J ] J)) yM • (6.3) 

n 

After a bit of algebra we obtain from (6.2) 

K ab (j) = 5 ab lim — ( V " ) =6ab dtpsvfcj) , (6.4) 

where p S v(£,]j) is defined as in (4.3) but with nonzero j. 

On the low-energy effective theory side we need to differentiate the log of the effective 
partition function 

Z% = J Vt L Vt R exp ( - j d 4 x £*£\ (6.5) 



with respect to j a and j b . Using (5.17) with mi ns t = and (5.19) we obtain 

$ 2 

ABCD 



K ab (j)= £ H tT(t a T A T B )tr(t b T c T D ) 



l - d'xd'y + + «&%)(y) 



: E Jftf ^(t a {T A , T B }) tr(t b {T A , T B }) 



AB "JAJB 



/\ conn 
d 4 x (vrf vrf (x)vrf vrf (0) + ^vrf (*)vr«(0)) j , (6.6) 

where "conn" denotes the connected part of the correlation function. To obtain (6.6) we 
have done the contractions A = C , B = D and A = D, B = C and symmetrized in A, B. 
The contraction A = B, C = D corresponds to disconnected diagrams and yields zero. 3 



37 The contribution of this contraction is proportional to Y^ b tv(t a T b T b ) oc tr(t a ) = since ^ 6 T 6 T b is 
the difference of the quadratic Casimir operators of SU(A r /) and Sp(iV/) and hence proportional to 1. 
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The dependence of (6.6) on j is contained in the masses of the NG modes. Evaluating the 
connected part in one-loop approximation 38 we find that it diverges for j — > 0, 



J d 4 *^£(*K^£(o) 



1 



d A p 1 
(2tt) 4 (p 2 +m\){p 2 +m 2 B ) ~ 16^ 



In 



A 



(6.7) 



where we have used (5.20) and A is the momentum cutoff of the integral, in which we 
have also absorbed /a,b and <3?h- Thus the dependence of the integral on A and B has 
disappeared, and (6.6) becomes 



$ 2 

K ab (j) ~ QabjTTz In 



A 



with 



Q<* = E JWW ^ a {T A , T B }) tr(t b {T A , T B }) 

AB 4 JAJB 



(6.9) 



where we have included a factor of 2 for the left- and right-handed NG modes in the loop. 
To evaluate Q ab we consider three cases: 



1. For A = B = 0, the contribution to Q ab vanishes since tr(t a ) = 0, 



Q 



(i) 

ab 



0. 



2. For A = 0, and A ^ 0, i? = 0, the contribution to Q a b reads 



^ab 



^ tr(t a T c ) tr(t b T c 



2<5qfe 



(6.10) 



3.11) 



where in the first equation the factor of 2 reflects the two possibilities above and in 
the second equation we have used tr(t a T c ) = S ac . 

3. For A and B ^ 0, the contribution to Q can be obtained from [84, eq. (48)] 
by replacing 2Nf —> Nf and multiplying by 8 to correct for the difference in the 
normalization of the generators (which in [84] is ti(t a t b ) = S a b/2). This yields 



off 



(N f -A)(N f + 2) 
2iV// 4 

Summing up these contributions, Q ab is given by 

'(N f -4)(N f + 2) 



(6.12) 



Q 



ab 



Qab + < 3ab' ) + Qab 



s, 



ab 



+ 



2iV// 4 ' iVj/o 2 / 2 

and our final result for the scalar susceptibility from the effective theory is 



(6.13) 



Kab(j) 



>ab 



{N f -4){N f + 2) 
2JV// 4 



+ 



H In 



16vr 2 



A 



(6.14) 



3 This is a valid approximation as we are interested in the infrared limit of the theory. 
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Now let us compare (6.4) and (6.14). First of all, note that the constant part p sv (0) 
in (6.4) does not contribute to K a t> since 

x £ 2 - i 2 
d t { p +j 2V = Q - ( 6 - 15 ) 







Hence only the difference p sv (£) — p sv (0) is relevant. To reproduce the singularity ~ ln(A/j) 
in (6.14), we must have 

PsAO - Psv(O) = for £>0 (6.16) 
in the vicinity of £ = 0, where C = p' sv (0). Now, 

d t (g2 +i 2 )2 ~Cln\ - 1 , (6.17) 
and thus the slope of the singular value density at the origin is given by 



o 



Psv(0) 



(N f -A)(N f + 2) + 2 



2N f f* N f f 2 p 



<J>2 

16vr 2 



(6.18) 



Note that <J>h, /o, and / are functions of p. To what extent this result is still valid at 
p < oo will be discussed in section 6.4. 

For Nf = 2 the Smilga-Stern method used above does not work since in (6.1) we then 
have J = jl so that K a b(j) cannot be defined as in (6.2). A similar problem occurs in 
the derivation of the slope of the Dirac eigenvalue density, where the Smilga-Stern method 
fails for Nf = 1. In that case the slope could still be computed using partially quenched 
perturbation theory [84], and it was found that the result obtained from the Smilga-Stern 
method remains valid for Nf = 1. It is therefore tempting to speculate that (6.18) remains 
valid for Nf = 2, but to confirm this we would have to compute p S v(£) in partially quenched 
perturbation theory. Such a rather complicated calculation is deferred to future work. 

Let us add two comments here. First, we could relax the assumption Jr = —Jl, and 
in particular we could set Jr (or Jl) to zero. 39 This would give us the slope of the density 
of the left-handed (or right-handed) singular values, which is 1/2 of the full slope because 
the factor of 2 mentioned after (6.9) would be absent. Second, we have now computed 
p sv (0) and Ps V (0) and therefore obtained information on the singular value density near 
zero. An analytical result can also be computed for asymptotically large £, which, owing 
to asymptotic freedom, can be described by the free theory without coupling to the gauge 
field. Thus the whole singular value spectrum can be understood at least qualitatively by 
an interpolation of two tractable limits. We obtain (for an arbitrary number A^ c > 2 of 
colors) 

^v(o^^e(e 2 + 2^ 2 ) for e^oo. (6.19) 

For p — > this reduces to twice the Dirac eigenvalue density in the free limit, as expected. 
An outline of the derivation is given in appendix E. 



39 This results in a projection on the topologically trivial sector (see (2.10)), which is immaterial in the 
p-regime. 
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6.2 Intermediate density 

The calculation at intermediate density is very similar to that at infinite density. Since the 
fundamental microscopic theory is unchanged, equations (6.1) through (6.4) also remain 
unchanged. On the effective theory side, (6.5) is replaced by 

Zl s = J VY> L VZ R VV exp (- J d 4 x Cl^j . (6.20) 

We could now go through a similar calculation as in section 6.1, using (5.24) through (5.26). 
However, it is easier to note that the only difference is the replacement of the two U(l) 
fields L and R by a single U(l) field V. The only contributions of the U(l) fields to Q ab 

(2) (2) 

are in Q ab , and it follows from the calculation in section 6.1 that we can obtain Q ab for 
the present case by dividing the result in (6.11) by 2. Everything else remains unchanged 
so that the slope for Nf > 4 is now given by 



Psv(0) 



(JV>-4)(iV) + 2) + 



2N f f* ' N f fgp 



$2 

(6 - 21) 



Again, $i, /o, and / are functions of \i. It is tempting to speculate that (6.21) remains 
valid for N f = 2. 

6.3 Zero density 

Let us now consider strictly zero chemical potential. Again, equations (6.1) through (6.4) 
remain unchanged, but the coset space of the effective theory is now S\J(2Nf)/Sp(2Nf). It 

(3) 

follows from the calculation in section 6.1 that Q a b is now entirely given by Q ab and that 

(3) 

we can obtain Q ab for the present case from the result in (6.12) by replacing Nf — > 2Nf, 
f — > F, and dividing by 2 since the left- and right-handed modes are already contained in 
SU(2JV»/Sp(2JV». This yields for N f > 2 

(N f -2){N f + l) *j_ 

P^>- NfF 4 167r 2' ^ ZZ ) 

where $l and F are now independent of ji. To what extent this result is still valid at 
nonzero [i will be discussed in the next subsection. 

6.4 Relation between the three results 

In the previous three subsections we have obtained three different results for the slope 
Ps V (0) for fi = oo, intermediate fi, and \x = 0, respectively. At first glance it does not seem 
possible to interpolate them smoothly, and we thus encounter a puzzle: How are the three 
results related? What actually happens to the spectrum if fj, is continuously changed? 
Below we argue that there is no puzzle here. To simplify the presentation we divide our 
discussion into three parts, the first one for three-color QCD at n = as an instructive 
model case for our problem, the second one for the low-density region, and the third one 
for the high-density region of two-color QCD. The arguments are analogous, however. 
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Figure 5. Left: Dirac eigenvalue density in three-color QCD at ji = for two massless flavors and 
one flavor with mass m, in units of toE 2 /327t 2 _F 4 . The intercept with the vertical axis is arbitrary 
due to renormalization [83]. Right: Slope of the density in units of £ 2 /327r 2 F 4 . The dotted lines 
in both plots correspond to the slopes {Nf = 2) and 5/3 {Nf = 3). 



6.4.1 Zero density {(5 = 2): A journey from Nf = 2 to 3 

So far we have argued that there are three effective theories for dense two-color QCD and 
that as a function of p they change smoothly from one to another. This situation has an 
exact counterpart in three-color QCD for Nf = 3 at p = 0, where the strange quark mass 
serves as a "knob" to interpolate between the chiral perturbation theories for Nf = 2 and 
Nf = 3. Therefore we study this simpler case first before considering the more exotic case 
of dense two-color QCD. 

The original Smilga-Stern relation [2], derived for the Dirac eigenvalue density (not 
the singular value density) in three-color QCD (/3 = 2) at /i = with Nf flavors in the 
chiral limit, reads 

S S 2 N 2 f - 4 

PW = -+ o 9 2F 4 M l A + °W ( 6 - 23 ) 
7r SLis^b^ Nf 

with the chiral condensate £ and the pion decay constant F. The coefficient of |A| depends 
on Nf. For example, the slope vanishes for Nf = 2 but is nonzero for Nf = 3. It is not 
clear from this expression alone how the slope changes if we add a nonzero strange quark 
mass to change the number of light flavors continuously. 

The generalization of (6.23) to nonzero degenerate masses was given in [82] 40 and later 
extended to nondegenerate sea quark masses in [83]. Therefore we can employ the results 
of [83] to study the density p{X) and the slope p'(X) for Nf = 3 with a massive strange 
quark. Setting (mi, m,2, ^3) = (0,0, m) in [83, Eq. (17)] we plot p(X) and p'(A) in figure 5 
as a function of X/m. 

The curves nicely interpolate between two limits: For A <C m the strange quark is heavy 
relative to the probed scale and we get the slope for Nf = 2: p'(X) oc (Nj — 4)/iV/ = 0. In 
the limit m ^ A (<C A ~ AirF) the strange quark is light relative to the probed scale and 

40 There is a typo in [82, Eq. (84)]. The second term in the square brackets of that equation must be 
multiplied by 7r, as is evident from [82, Eq. (83)]. 
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we get the slope for Nf = 3: p'(A) oc (iV| — 4)/iV/ = 5/3. The transition occurs smoothly 
around A ~ m. Thus we can draw the conclusion that no contradiction arises from different 
values of the slopes for Nf = 2 and Nf = 3, because at nonzero m they correspond to 
different domains of the spectrum. 

This finding can be interpreted within partially quenched chiral perturbation theory 
as follows. In this method we add valence flavors and compute the spectral density from 
the valence quark mass dependence of the chiral condensate [85] . We therefore deal with 
two classes of mesons, one being made of only sea quarks, and the other being made of 
valence quarks (and sea quarks). If the latter ("valence mesons") are much heavier than 
the former ("sea mesons"), i.e., A>m, then all three sea flavors contribute to the valence 
quark mass dependence of the chiral condensate, implying Nf = 3. Conversely, if A <C m, 
the heavier sea mesons are decoupled, reducing the computation to Nf = 2. The transition 
between these two cases occurs around A ~ m, i.e., when the masses of the sea and valence 
mesons are roughly equal. 

6.4.2 Low density 

In this subsection we discuss the relation between the results (6.21) and (6.22) for |i< A. 
The chemical potential plays exactly the same role as the strange quark mass in the previous 
subsection, both acting as explicit symmetry-breaking parameters. This analogy is the basis 
of our following argument. 

We first note that p S v(£) can be computed in partially quenched chiral perturbation 
theory, starting from 



and setting j v = j' v = i£ + e (with e — > + ) at the end of the calculation. We will not 
actually perform this computation but use (6.24) for qualitative estimates. In comparison 
to the usual setting [82, 84], j and j v ,j' v ~ £, correspond to the sea and valence quark 
masses, respectively. In the following j is always assumed to be infinitesimal. Our main 
concern is the competition between /x and j v . 

We can now replace Z Nf+2 \2(j', jv, j' v ) by an effective partition function formulated in 
terms of NG modes. We have two options, either the partially quenched extension of theory 
L, or that of theory I. Let us discuss them separately. 

• Theory I: For the partially quenched extension of theory I to be valid, the condition 
(i) of section 5.5 must be satisfied for all NG modes, i.e., their masses must be much 
smaller than rri£. The new ingredient here is that, in addition to the NG modes 
discussed in section 5.4, we now also have NG modes containing the valence quarks 
corresponding to j v and j' v , which we call valence NG modes. At low density their 
masses are of order \/j v A ~ v^A, see (5.10a) with F ~ <I> L ~ A, and hence the 
condition is 




(6.24) 




(6.25) 
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In this domain the slope p' sv (0 is given by (6.21) at leading order of the low-energy 
expansion. 

• Theory L: For the partially quenched extension of theory L to be valid, the masses 
of all NG modes must again be much smaller than mg. This time the masses of the 
valence NG modes are of order y/j v A ~ \/£A and \J j v A + p 2 ~ see (5.10) 

1/3 

with F ~ <I> L ~ A. The condition is therefore 

£<A. (6.26) 

This is not end of the story, however: Theory L is more complicated than theory 
I, because we have two scales p, and \/j v A (analogous to m and A in section 6.4.1) 
whose ratio controls the final result in a nontrivial way. Based on our experience in 
section 6.4.1 we expect the following. 

— For \i <C y/ j v A ~ \/£A : All sea NG modes contribute, and the result for p' sv (^) 
agrees with the result (6.22) at /i = 

— For fi 3> \/j v A ~ y/£A : The NG modes with masses of order \x decouple from 
the computation of p S v(£)> anc ^ theory L reduces to theory I in which the heavy 
modes have been integrated out. The slope thus agrees with (6.21) from theory I. 

Putting everything together, we see that in the regime p <C A the results from theory 
I and L for the slope p' av (£) are valid in the following domains: 

u 2 

(6.21) from I : f < j- , (6.27) 

u 2 

(6.22) from L : fj- < £ < A . (6.28) 

These findings are illustrated in figure 6 (left). The slope is first given by (6.21), and we 
conjecture that it changes smoothly to the value given by (6.22). To avoid confusion we 
point out that for p = the window in which the slope is given by (6.21) shrinks to zero 
so that the slope at the origin is given by (6.22). 

We now present a nontrivial cross-check of our conclusion that does not hinge on 
partially quenched chiral perturbation theory. Let us return to our discussions in sections 
6.1-6.3 without valence quarks and consider taking the limits fj, — > and j — > while 
keeping the condition p? S> jA. In this case the two types of NG modes have masses 
of order \fJK (lighter) and \x (heavier), respectively. The one-loop integral (6.7) in the 
effective theory can be evaluated for the cases of two lighter, two heavier, or one lighter 
and one heavier NG modes circulating around the loop. For these cases we obtain 

A I A 

d P To Tv> ~ ln ~ . ( 6 - 29 ) 

y (p 2 +jA) 2 j 

A 1 A 2 

d 4 p t^, o\~o ~ In — o ; (6.30) 

{p 2 + p 2 ) 2 p 2 K ' 
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Psv(£) * 




slope from (6.22) (L) 



slope from (6.21) (I) 
1 



slope from (6.18) (H) 



\ 



slope from (6.21) (I) 



P 2 /A 



gm 2 ,/A 



Figure 6. The behavior of the singular value density near zero. Left: low density (fi <C A), right: 
high density (so that m,y -C A), see text for details. 



d 4 p 



fi 2 In 4-- jA la- 



in ■ 



A 2 



(p 2 + jA)(p 2 + /i 2 ) [i 2 — jA fi 

Therefore the sum of all one-loop contributions to is given by 

i A oi A2 
a In h p m 

j P 



for /j 2 > j A . 



(6.31) 



(6.32) 



with prefactors a and f3 that also contain traces of the generators. The two infrared 
singularities in (6.32) (generated by j — > and fj, — > 0, respectively) must be matched by 
corresponding singularities in the microscopic theory. Motivated by figure 6, let us assume 
that p S v(0 can be approximated by a straight line with slope a for £ < £ c and by another 
straight line with slope a + (3 for £ > £ c , with £ c an unknown function of /x. Then (6.4) 
becomes 



d£p a v{£) 



2\2 



d£ (a£ + const.) 



+ 



d£ [(a + /3)£ + const.] 



(e+j 2 ) 2 



, f c . , A , A _ , A 
~ a In h (a + p ) In — = a In h p in — . 



Matching (6.32) and (6.33) yields 



A 



(6.33) 



(6.34) 



in agreement with the argument based on partially quenched chiral perturbation theory. 
6.4.3 High density 

We now clarify the relation between (6.18) from theory H and (6.21) from theory I at large 
fi. The arguments are analogous to those at small /j,, except that mi ns t now plays the role of 
[j, as an external symmetry-breaking parameter. First of all we require that the masses of 
all NG modes (sea and valence) must be sufficiently below rri£. At high density the masses 
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of the valence NG modes in the partially quenched theory are of order \J j v Aj g ~ y^A/g 



and ^j v A/g + mf nst ~ ^A/g + m^., see (5.20) with / A ~ p [74] and $ H ~ p 2 &/g 
[81], where 5 is the running coupling constant. Using the relevant results for mi in table 1 
we obtain the following bounds on the values of £ below which p S v(£) can be computed 
from the partially quenched extensions of the effective theories I or H in the regime where 
m v > <C A: 

ami, 

£ < for theory I , (6.35) 

£ < gA for theory H . (6.36) 

Thus Ps V (£) is given by (6.21) from theory I in the range (6.35). We note that the scale 
grn?,/A goes to zero rapidly as p — > 00. 

On the other hand, theory H has two scales, \J j v Aj g and mj ns t. A rerun of the 
arguments at small p then shows that in the regime where m„/ <C A the slope p' sv (£) is 
given by the results (6.18) and (6.21) in the following domains, 

gml, 

(6.21) from I : £ < —3- , (6.37) 

gml, 

(6.18) from H : —^L < £ < 5 A . (6.38) 

This is illustrated in figure 6 (right). For fx — > 00 the window in which the slope is given 
by (6.21) shrinks to zero so that the slope at the origin is given by (6.18). 

The second argument presented in section 6.4.2 works in exactly the same way here 
and leads to 

gml, 

tc~-^~ (6-39) 

as expected. 

Note that the result (6.21) is actually valid for all < fi < 00. We can replace p' sv (0) 
by Ps V (£) in (6.21) for sufficiently small £, the upper bound of which is given in (6.34) 
and (6.39) at small and large p, respectively. At intermediate density we do not have an 
estimate for the upper bound because me is unknown in this region. 

So far we have explained what we believe is the most reasonable behavior of the 
singular value density at nonzero p based on the partial quenching technique and the 
analogy to three-color QCD at p = with a heavy strange quark. For a solid proof of our 
conjecture shown in figure 6 one would have to compute the slope in partially quenched 
chiral perturbation theory explicitly, but this is beyond the scope of this paper. It would 
be interesting to check our new Smilga-Stern-type relations by lattice simulations. This is 
possible in principle as the infamous sign problem is absent in this theory. 
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7 Finite-volume analysis: Leutwyler-Smilga-type sum rules and random 
matrix theories (f3 — 1) 

7.1 The e-regime 

In this section we study two-color QCD with diquark sources in a finite volume V4 = L 4 
(and again in the chiral limit). As in QCD there is a regime, the so-called e-regime [3], in 
which the kinetic terms in the effective chiral Lagrangian can be neglected so that the theory 
becomes zero-dimensional and the partition function is dominated by the zero-momentum 
modes of the NG particles. The condition for the e-regime is 

— «L<^— , (7.1) 
me m NG 

where me is again the mass scale of the lightest non-NG particle and 7Bng is the mass 
scale of the NG particles that are included in the effective theory. The first inequality 
in (7.1) means that the contribution of the non-NG particles to the partition function 
can be neglected, while the second inequality means that the Compton wavelength of the 
NG particles is larger than the size of the box, which in turn implies that the functional 
integral over the NG fields can be replaced by a zero-mode integral over the coset space 
parametrized by them. For the three different effective theories in section 5 the values of 
me are given in section 5.5, and the mNG are given in (5.10), (5.20), and (5.27). Note 
that for zero diquark sources and finite L the second inequality in (7.1) is always satisfied 
in theory I since tting = 0. Note also that while the domains of validity of the effective 
theories overlap (see section 5.5), this is not the case for the corresponding e-regimes: 
At low density it follows from (7.1) that l//x <C L for I and L< 1//J for L, and these 
conditions are mutually exclusive. Similarly, at high density we have l/m^i <C L for I and 
L <C l/mrf for H. 

At [i = 0, it is well known that there exists a scale Et below which the eigenvalue 
spectrum of the Dirac operator obeys chiral random matrix theory [82, 86]. The scale 
Et is called Thouless energy, borrowing the nomenclature from mesoscopic physics. The 
equivalence between random matrix theory and the zero-momentum limit of the partially 
quenched effective theory shows that Et can be understood as the energy above which 
the condition (7.1) no longer holds in the partially quenched theory and the modes with 
nonzero momentum start to contribute. 

We now comment on the Thouless energy for the singular value spectrum at < fx < 00, 
based on the partially quenched extension of the three effective theories introduced in 
section 6.4.2. For simplicity we let j = in (6.24) and concentrate on the "spectral mass" 
j v . We assume that the condition L S> I /me is satisfied for all cases considered below. 

• Theory I (0 < \i < 00): From (5.27) we find that the masses of the valence NG modes 
are given by m> v -^Q ~ j v 

^i// 2 ~ C&i/ f 2 , where all low-energy constants depend on \i 
implicitly. 41 For m v NG <C l/^ the j„-dependence of the partition function is governed 

41 We are sloppy about the distinction between / and fo here, but they are of the same order of magnitude 
so that the distinction does not change our discussion. 
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by the NG modes with zero momentum. Thus the Thouless energy is determined by 
VEt^i/P = \ -> E t = J^ i - (7-2) 

• Theory H {mrf <C A): From (5.20) we have m^ NG ~ j v ^u/f 2 ~ (see section 
6.4.3) and j v ®n/f 2 + mf nst ~ j v A/g + m? nst . Assuming m inst « 1/L we obtain 

^K/g=~ -> E T = ^- (7-3) 

• Theory L (^x < A): From (5.10) we have m- 2 NG ~ j^l/^ 2 and jt^iV-F 2 + (2^) 2 , 
where this time the low-energy constants are //-independent. Assuming fi <C 1/L we 
obtain 

V^l/F 2 = i -> £ r = _^. (7.4) 

Thus in all cases Et oc 1/v^4- Essentially, this is due to the fact that the diquark 
source enters the effective Lagrangian linearly at any [i. This completes our discussion of 
the Thouless energy. 

In the e-regime we can compute exact sum rules for the inverse singular values of the 
Dirac operator. This will be done in the following three subsections for the three different 
density regimes. For technical reasons we start with intermediate density this time. We will 
also derive chiral random matrix theories that allow us to compute microscopic correlation 
functions of the singular values. 42 The sum rules are simply moments of these correlation 
functions. The results for p sv (£) from the random matrix theories are valid in the range 
i < E T . 

7.2 Intermediate density 

In analogy to the analysis of Leutwyler and Smilga [4] we first project the partition function 
onto sectors of fixed topological charge and then expand it in powers of the diquark sources. 
This is done both in the microscopic theory and in the effective theory. Matching the 
coefficients of the sources then yields sum rules for the inverse singular values. 

On the QCD side, we start from (2.10) and assume that there are no accidental zero 
modes. Without loss of generality we therefore set til = and ur = v > 0. 43 We also 
introduce the notation 

(Odet'ptz^/ 2W 
{Uh (det'pt^^W ' 1 J 

where the subscript v indicates that the average is only over gauge fields with fixed topo- 
logical charge v and the prime, as always, means that the zero modes are omitted. Using 

det' {rfD + J{J L P L + J r JrPr) = [J' det(£L + + JrJr) (7-6) 



42 In addition, a chiral random matrix theory for QCD with isospin chemical potential (/3 = 2) will be 
derived in appendix F. 

43 As usual, for v < the final results of this section remain valid, except that v must be replaced by \v\. 
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the partition function for fixed topology is given by 

ZAJu Jr) - [pf(4)Pf(j fl r/n'det i/2 U+ JrJr 



Using the formula 



(7.7) 



det(l + e) = 1 + tr e + 1 [(tr ef - tr(e 2 )] + 0(e 3 ) (7.8) 



2 

we expand Z„ in powers of the diquark sources, 
Z u (Jl, Jr) 



(det'( J Dt jD) ^/2) YM ^p f(J t )Pf(J/? 
= 1 + 



tr JIJl/y>±\ tr ( J l J L? (yf 1 \ | (trjjj L ) 2 //^/ 1 



2 \ ^ ^ / 4 \ ^ f A / ' 8 \ \ ^ £ 2 



+ (L o R) 



tiJ ] L J L txJ ] R J R /(^ 1 1 



On the low-energy effective theory side, let us first consider the case Nf > 4. We start 
from (5.24) and neglect the kinetic terms to obtain the finite-volume partition function 

Z eS (J L ,J R ) = J dZ L dZ R dV exp [U 4 $iRetr(J L £ L - J R E R )V] , (7.10) 

where the integrals over Ex R and V are no longer functional integrals but simple integrals 
over SU(A r j)/Sp(A r j) and U(l), respectively. It is convenient to define 

J i = J i V^ l (i = L,R) . (7.11) 

As explained after (2.12), a nonzero 6>-angle can be introduced by redefining, e.g., Jl — > 
j L Q~ ie / N f anc l J R — y J R e ld / N f. To project onto topological sectors we note that the 
inversion of (2.12) is 

Z v = r d9 e- tud Z(6), (7.12) 
^ Jo 

which in the present case gives 

dZ L dZ R dV — exp [ - iv0 + Retr( J L e~ ie / N f £ L - J R e id / N f E R )V] 

/df) 
dZ L dZ R dV — exp [ - iNfvO + Re tr(J L e~ ie T, L -J R e ie E R )V] 

where in the second step we have redefined — > Nf0 and used the periodicity of the 0- 
integral to extend the integration region from [0,2w/Nf] back to [0, 2n]. In the last step 
we have introduced the new U(l) integration variables 



L = e~ t0 V and R = e ie V . (7.14) 



_ 4i _ 



Parametrizing Ej = UiWf with U { G SU(JV»/Sp(JV» as in (5.11) we note that 

Retr(J R i?E fi ) = Retr(jJi? t S^) = Retr(- J R R?U R IU R ) . (7.15) 

Redefining the integration variables R FV and — >■ C/^ and combining L and i? with 
Si and £/j as in (5.18) we thus obtain from (7.13) 

Zf(J L ,J R ) = J dELcERdLdRiLR)^" exp [Retr(J L L£ i + J R RE R )] 

= J dt L dt R det^ 2 (t L t R )exp [Retr(J L E L + J R t R )] . (7.16) 

Recalling Ej = UiIUj with C/j S U(A r j)/Sp(A r j) we note that we can extend the integration 
to be over Ui E U(iVy) because the generators of Sp(iVj) leave / invariant and therefore 
drop out of the combination UilUj '. We now define, for a complex antisymmetric matrix 
X of even dimension Nf, the function [87] 

9v (X) = [ dU(detU) u exp[Retr(X^UIU T )], (7.17) 

Jv(N f ) 

in terms of which we have 

Zf(J L ,J R ) = g u (Jl)g u (J R ). (7.18) 

The expansion of the integral (7.17) in powers of X was computed in [87]. Adapting eq. (40) 
of that reference to our case, we have with a = Nf + v — 3 

g v {X) =N{Y>iX) u \l + A a tTX ] X + B a {tiX ] X) 2 -C a ti{X ] X) 2 + 0{X Q )\ (7.19) 
with coefficients 

' B <* = -o 7 "t} , , . C a = A , ■ * ■ - (7.20) 



2(a + 2)' 8a(a + 2)(a + 3) ' 4a(a + 2)(a + 3) 

and a normalization factor AA that depends on Afy and v but is not important for our 
present purposes. (A" can be determined from [88] if desired, and for v = we have A" = 1 
for all Nf.) Equations (7.18) and (7.19) imply that in the limit Jl/ r — > only the v = 
sector survives. This is analogous to QCD at fj, = 0, where in the chiral limit the topological 
susceptibility vanishes [89]. We thus obtain 

Z '" U{ Jl ' J/>> ' " 1 + [A a tr j[ J L + fl a (tr 4Jl) 2 - C a tr( J| J L ) 2 + (L o fl)] 

+ A 2 (tr f L J L )(tT J R J R ) + 0{ J 6 ) . (7.21) 



AA2[Pf(jt) P f(J H )]^ 



We can now match the right-hand sides of (7.9) and (7.21) to obtain sum rules for the 
inverse singular values, 

E'^) = <E'^> ((E^)) -((E^) 2 ) 

(7.22a) 
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/2 C4 



n ^Rn I v 



y — 

m ^Rm 



y — 

Z^ C2 

n ^Ln 



(7.22b) 



where A a = {V^i) 2 A a , B a = (V4<3?i) 4 -B a , and C a = {V^i) A C a . These sum rules imply 

2\ 



Note that the conditions 



Z_> C2 



E 



/ 1 



> 



E 



/ 1 

t 2 



and 



+ 165 



Ei 1 1 
C2 Z-> f2~ 

n ''Ln n ^Rn 



E'^> = 8 ^- ( 7 - 23 ) 



16£ a - 8^ > o 



(7.24) 



both lead to the inequality 2-B Q > A^, which is satisfied nontrivially for all Nf > 4 since 



2B a - Ai 



2a(a + 2) 2 (a + 3) 



> 0. 



(7.25) 



The second condition in (7.24) also makes it clear that in general the left- and right-handed 
sums are different (for a fixed gauge field). Higher-order sum rules can be computed by 
expanding the partition functions to higher order in the diquark sources. 

Let us now consider Nf = 2 and for simplicity take j R and jl real. In that case (7.9) 
becomes 



(-jur) 



Ln 



Zw £2 



Rrn / I v 



+ 



y^ c2 



+ 3i 



y — 

Z^ t2 c 



+ (L^R) 



+ 0(f 



y m<n ^Lm^Ln I v 

The finite- volume partition function obtained from (5.28) in the e-regime is 



(7.26) 



Z cS (j L , or) = [ dV exp [ - 2 (j L - j R ) Re V] 
Jv(i) 

with ji = jiV^i. Introducing a #-angle as before we thus have 

Zf(jLjR) = J dV £ exp [ - ivQ - 2 Re(j L e"^ 2 -j R e w ' 2 )V] 

= [ dLL u exp[-2j L ReL] I dR (R r ) u exp [2j R Re R] 
Ju(i) Jum 



(7.27) 



'TJ(l) 
Iv{-2j L )U2j R ) 



(7.28) 



in analogy with the steps that led to (7.13). In the last line we have recognized the integral 
representation of the modified Bessel function l v . Expanding in the diquark sources gives 



Zf(j L ,jR) 
{^■)- 2 (-jL3RY 



y _|_ ji j_ jl _|_ 3l3r _|_ 



Jl+Jr 



u + 1 (v + 1) 2 2(i/ + 1)(i/ + 2) 



+ 0(j e 



(7.29) 
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and matching (7.26) and (7.29) yields the sum rules 

r'^\=/v'l\= //V' 1 Vr ; M\ = (F4#l)4 

(7.30a) 

E'^V) -(E'jrr) < 7 - 30b » 

m<n ^Rm^Rn/u \ m<n ^Lm^Ln / u ^tlj^t^ 

Note that the sum rules in (7.30a) follow from those in (7.22) by setting iVy = 2. 
We observe that there is a "decoupling rule" for Nf > 2, 

E^)(£^)H£^)(£^). <-) 

m >i?m / \ n ^Ln /IV \ n ^Rn I v \ n ^Ln I v 

A similar decoupling of the left- and right-handed modes is expected to hold for higher 
moments as well since the finite- volume partition function factorizes, see (7.18). It is 
therefore sufficient to consider only one of the factors. Considering Nf > 4 for the time 
being, we divide (7.7) and (7.18) by [Pf j\] v and then take the limit Jl — > to obtain 

dU (det U) v exp [V&i Reti:(J R UIU T )} oc [Pf J R ] U (TT'det 1/2 (1 + j \ 

'V(Nf) " \ n V Z Rn ) j v 

(7.32) 

modulo an unimportant normalization factor. This has the same form as the mass- 
dependent finite- volume partition function at [i = for Nf > 2 [88], 



Zf(M)= [ dU(detU) u exp 

Jv(2N f ) 



-V 4 ZRetr{M j UIU T ) 



«(P<Mr(n^ /2 (* + ^-)) M-[_ U "), (7.33) 



\A„>0 




where M is the quark mass matrix, E is the absolute value of the chiral condensate, the 
i\ n are the Dirac eigenvalues at \i = 0, and the normalization factor was omitted again. 
This implies the correspondence 

{$i, J R , N f , t Rn }^ <— > {|E, M, 2N f , X n (> 0)} m=0 (7.34) 

and similarly for R -R- L. For Nf = 2 this correspondence remains valid since the U(2)- 
integral in (7.33) then gives I v {_— mV^E) with the quark mass m, which is to be compared 
with I v (2j R Vt$i) in (7.28). 

It is well known that to lowest order in the e-regime the system can alternatively by 
described by chiral random matrix theory. For our present case, i.e., two-color QCD at 
intermediate density in the chiral limit and in the presence of diquark sources, we obtain 
the chiral random matrix theory 44 

zr T ( Jl, Jr) = 9 T T (JI)9T T (Jr) (7-35) 



4 Random matrix theories for singular values are also considered in a very different context in [90, 91]. 
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with 



RMT 



(J) = fdA e -"*<^> Pf J ® lN T , lN < ® A 
J \t Nf ®{-A T ) J®t N+u J 

N f /2 

= fdA e~ ff ^ J] j? det (AA T + |j/| 2 ) , (7.36) 
J f=i 

where A is a real matrix of dimension N x (iV + u). In the second line, we have applied an 
orthogonal transformation O to bring the antisymmetric matrix J to the standard form 

J = o{ A ° . diagGi,...,^)^^ 
\-diag(ji, . . . ,j Nf / 2 ) y 

where the jf can be complex. For u < we have to replace by (j*)'"' and ^4^4^ by A T A 
in (7.36). Note that the second line of (7.36) clearly exhibits the correspondence of the 
random matrix theory for our present case and for two-color QCD at \i = with nonzero 
masses and without diquark sources [51]. 

Going through the standard steps of converting the random matrix theory to a sigma- 
model we find in the limit N — > oo that g^ MT {J) = g u (V2NJ), which shows that the 
random-matrix partition function is equivalent to the finite-volume partition function, 45 
provided that the dimensionless random-matrix diquark sources and singular values (i.e., 
the square roots of the eigenvalues of A T A) are related to the physical quantities by 

J i = JiV^i/V^N and ii = iiV^i/y/2N (i =L,R). (7.38) 

Note that the partition function only factorizes in the chiral limit. If we include quark 
masses the random-matrix partition function is 

/ J L A L -M T \ 



Zf MT ( JL, Jr, M) = J dA L dA R e -Ntr(AlA L +A%A R ) pf 



-A\ j{ -Mt 
M -,P R -A R 
\ M* A T R -J R ) 



(7.39) 



where Aj_, and Ar are again real N x (N + v) matrices. This is a natural extension of the 
random matrix theory constructed at high density in the absence of diquark sources [16, 93]. 
How the dimensionless random-matrix quark masses are related to the dimensionful masses 
depends on the density. At high density we have M = My / 3T4/./VA/27r [16]. What sets 
the scale for the masses at intermediate density is a dynamical question that cannot be 
answered with the methods we employ here. 

For jj, = 0, the microscopic spectral correlations of the Dirac eigenvalues have been 
computed in chiral random matrix theory [94, 95]. The microscopic scale is defined by 
z = AV4S, and as an example we quote the microscopic spectral density [96] 

p?""(z) = Z 2 [J 2 a {z) - J a+1 (z)J a ^(z)] + \j a {z)[l - J Z dwJ a (w)] , (7.40) 



45 This equivalence is expected to extend to the microscopic correlations of the singular values, although 
this still needs to be proven using the partially quenched theory, similar to [92]. 
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where a = 2Nf + \i>\ and J denotes the Bessel function. According to the correspondence 
(7.34), the microscopic density of the Dirac singular values for a given chirality, 

Psv f, "(x)= lim (S^ 5(x - x n )) with x n = £ Rn V 4 $i or £ L nV±<$>i , (7.41) 

n 

is given by 

PsJ' u (x) = 2pl Nf ' u (2x) (7.42) 

for Nf > 2, and analogously for higher-order correlation functions. By construction all 
sum rules derived in this section are moments of suitable microscopic correlation functions. 
We have checked numerically that this is indeed true for the moments of the microscopic 
density of the singular values. 

Note that for Nf = 2 we could not obtain sum rules for, e.g., \Yl n ^-/^ Rn ) v or 
(X^n ^-/^Ln) v f rom t ne expansion of the partition function. However, we can obtain them 
as moments of the microscopic density. It follows from the Taylor expansion of (7.40) that 
the above two sums diverge for v = and 1, which is consistent with the observation that 
formally setting Nf = 2 in (7.22b) gives a negative or infinite result for v = or 1, respec- 
tively. But for v > 2 these sums converge, and the result is identical to the one in (7.22b) 
with Nf = 2. Since in the random matrix theory there is nothing special about Nf = 2 we 
expect that all sum rules derived for Nf > 4 remain valid for Nf = 2, if convergent. 

7.3 High density 

At asymptotically high density topology is suppressed and mi ns t = 0. In that case we 
expect to obtain the v = subset of the results of the preceding subsection. However, 
as long as mi nst is still nonzero the v ^ sectors are not completely suppressed. These 
expectations will be confirmed below. 

On the QCD side nothing changes, i.e., we match to (7.9). On the low-energy effective 
theory side, let us again start with Nf > 4. Neglecting the kinetic terms in (5.17) gives 

Z eS (J L ,J R ) = — ^- / dL L dE R dLdR exp [Retr(J L LS L - JrR^r) + KRe(tf R) N f /2 } 
Iq{k) J 

(7.43) 

with Ji = JiV^ n and k = 2Vj$m 2 ms jNf > 0. 46 The integrals over L,R are over U(l), 
while those over Sx,^ are over SU(A r j)/Sp(A r j). The normalization factor 1/Iq(k) has been 
added to ensure Z e ff(0,0) = 1 as in (7.10). Introducing a #-angle as before and projecting 
onto topological sectors using (7.12) we obtain 

Zf(J L ,JR> = [ dY, L dY, R dLdR^- 

x exp [ - ivQ + Retr( J L e~ id / N f LS L - J R e w / N f RZ R ) + KRe(L^R) N f/ 2 ] 
= JiE L iE R dL dR (L flt) a"/- exp [Retr(J L LE L - J R RZ R>] 

46 Although fo ~/iwe have k — > for fi — > oo since mf nat goes to zero much faster than 1/fJ 2 [67]. 
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/o(k) 



/d0 
— exp ( — iv9 + k cos 9) 



[zf(J L , J R ) of theory i] x |M , ( 7 .44) 



where in the second line we have first redefined L — > Le t9 / N f, R — > Re~ t9 / N f and then 
e —2i6/Nf _^ e -2i8/N f _^^jt ; anc [ i n the last line we have compared with (7.13). Prom (7.44) 
we can draw several conclusions. First, all Leutwyler-Smilga-type sum rules in theory 
H are identical to those of theory I (except for the replacement $i —> $h) since the 
relative factor I u (k) / Iq(k) is independent of Jl/r and therefore drops out when computing 
expectation values in sectors of fixed topology. Second, in the limit k — > oo the relative 
factor goes to 1 for any v , while for finite k the sectors with v ^ are suppressed and 
disappear completely for k — > as expected. Third, even in the presence of the anomaly 
the finite-volume partition function still factorizes into a left- and a right-handed part as 
in (7.18). 

For Nf = 2 the argument goes through in exactly the same way. 

Since the sum rules for theories I and H are identical, it is natural to expect that the 
random matrix theory and the microscopic correlation functions of the singular values for 
H are identical to those of I, with <I>i -R- <1>h (again, this would have to be proven using 
the partially quenched theory). The only difference is in the summation over topological 
charge. If, e.g., one wants to sum the microscopic spectral density over all sectors, one 
needs to take into account the ^-dependent relative factor in (7.44), i.e., 

£ %v(x) Zl £ %v(x) Zl /„(«) 

pl(x,6)= v ——— and fefl)= " . - — , (7.45) 

V V 

where we have suppressed the arguments Jl and Jr in Z\ and p sv . Note that in the limit 
of zero diquark sources Z\ = for v ^ so that p sv (x, 9) = Pg V (x) for both I and H. 

7.4 Low density 

On the low-energy effective theory side, we use (5.4) and for simplicity take Jr = — Jl = jl 
with real j. Neglecting the kinetic terms we obtain the finite- volume partition function 

Z cS (p,j) = J cffi exp [^ 2 FV 4 tr(S J BS t S) + jV 4 $i Retr(e~ ie/N f S d S)] , (7.46) 

where the integration is over SU(2iVj) /Sp(2Nf), the £> 2 -term has been absorbed in the 
normalization, and a #-angle has been introduced by J R -»■ Jr e i0 / N f and J L -> J L e~ ie / N f . 
The expansion of (7.46) in powers of j for arbitrary Nf is formally possible, but the 
analytical calculation of the expansion coefficients is a challenging mathematical problem 
which we do not address here. Instead, we will obtain partial results for the special case 
Nf = 2. Because of the group isomorphisms SU(4) ~ SO(6) and Sp(4) ~ SO(5) we 
can regard the coset as SO (6) /SO (5) ~ S 5 [88]. This approach has been explored in 
detail [71], and in the following we use the formulation of that reference. Adapting the 
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unnumbered equation just before eq. (7) of [71] to our case, we have modulo a j-independent 
normalization factor 



Z eS (v,j) = [ dnexp 



z(n\ + n|) + 2wn 4 cos 



with 



(7.47) 



z = 8/i 2 F 2 V4 and w = 2jV^ L 
Projecting onto fixed topology using (7.12) we obtain 



(7.48) 



s-< 



d9 e 



-iv9-\-2wri4 cos -, 



U(l) 



/ dn e ^ n l +n D I 2l/ (2wn 4 ) 
Js 5 



s s (2i/)! 



2// 



i + {um ^ + 0(f) 



(7.49) 



As before, we always assume v > 0. For v < we have to replace v — > \ v\. 
On the QCD side we again use (7.9), which now becomes 



Z u (u,j) oc j 



■2v 



1 



+ 0(f 



(7.50) 



where we omitted a j-independent normalization factor. Comparing the ratio of the j 2u 
and j 2u+2 terms in (7.49) and (7.50) we obtain the sum rule 



y- 

„ Sn 



2u + l 



2 J S 



dn e 2 ( n 2+ n D nf +2 



dfi e 2 ( n 2+ n D nf 



(7.51) 



For z = 0, (7.51) reproduces the sum rule in [88], as it should. 47 In the opposite limit 
z — > oo, the r.h.s. converges to 2(V4 ( I ) l) 2 / (^+1) in agreement with (7.30a) (with <I>i <-> 
The analytical calculation of the r.h.s. for finite z is nontrivial, and we only consider the 
special case v = 0, for which we obtain 



y- 



2(F 4 $l) 



-z-1 



+ 1 



(7.52) 



the plot of which is shown as a function of z in figure 7. The proof of this result is given in 
appendix G. The sum rules for the right- and left-handed singular values are also given by 
(7.52) but with the r.h.s. divided by 2. Higher-order sum rules can be obtained as usual 
by expanding (7.49) and (7.50) to higher orders in j. 



47 In this comparison we need to be careful. First, our $l corresponds to E/2 in [88]. Second, all singular 
values are doubly degenerate (for /3 = 1) at fx = 0, which means that the sum rule in [88] should be 
compared to our result divided by 2. 
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Figure 7. /x-dependence of the spectral sum rule (7.52) in theory L in the sector v — 0. The curve 
converges to 1/3 as fi — > and to 1 as fi oo. 



As before, to lowest order in the e-regime we can also describe the system by a random 
matrix theory, which at low density we can obtain by adding diquark sources to the known 
two-matrix model [18, 97], resulting in 

Z^(H,J L ,J R ,M) = 

( J L C-fiD -M T \ 

C T + fiD T j{ -Mt 



/ 



dCdD e-«"<(WD)pf 



M -J R -C-fiD 

\ M* C T + fiD T -J R ) 



(7.53) 



For earlier approaches, see [98, 99]. In (7.53), C and D are again real iV x (N + v) matrices. 
Converting (7.53) to a sigma-model we find in the limit N — > oo that the random matrix 
parameters are related to physical quantities by 

fi 2 = 2fi 2 F 2 V 4 /N , M = MV^ L /2N, J i = J i V A ^/2N, & = &y 4 $ L /2iV , (7.54) 

where in the chiral limit the £l (£r) are the singular values of C — fiD (C + fiD). Note 
that for fi ^ 1 the partition function does not factorize even in the chiral limit. However, 
for /t = 1 ("maximum non-Hermiticity" ) (7.53) reduces to (7.39), and if the chiral limit 
is taken it factorizes again. The computation of the microscopic correlation functions of 
the eigenvalues and/or singular values from this random matrix theory is a complicated 
mathematical task which we do not attempt here. 

In section 7.1 we have seen that the e- regimes of the effective theories L and I do not 
overlap. Nevertheless, as just noted, the random matrix theory (7.39) for I is the ft — > 1 
limit of the random matrix theory (7.53) for L. This is consistent with the observation that 
the sum of the sum rules for l/£ Rn and l/££ n in (7.30a) is the z — > oo limit of (7.51) for 
all v (with $i f-> ^l)- We are therefore tempted to conjecture that all sum rules of I are 
the z — > oo limits of the sum rules for L. 
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8 Conclusions and outlook 



In this paper, we have studied the singular values of the Dirac operator in QCD-like theories 
with Dyson indices /3 = 1, 2, and 4 at nonzero chemical potential [i. We pointed out that 
the Dirac singular values are always real and nonnegative and that the scale of the singular 
value spectrum is set by the diquark condensate at any \x. This is in contrast to the Dirac 
eigenvalues, which spread into the complex plane at nonzero \i and whose scale is set by 
the chiral condensate at small \x [6] and by the BCS gap at large \i [15, 16, 18]. We 
derived Banks-Casher-type relations for all three values of (3 and then concentrated on the 
(3 = 1 case, for which we identified three different low-energy effective theories with diquark 
sources at low, intermediate, and high density within the whole BEC-BCS crossover region, 
and clarified how they are related to each other from the point of view of integrating out 
heavy degrees of freedom. We derived exact results, such as Smilga-Stern-type relations and 
Leutwyler-Smilga-type sum rules, which (together with the Banks-Casher-type relation) 
concern the connection between the singular value spectrum and diquark condensation. We 
have also identified the e-regimes of the effective theories and constructed the corresponding 
chiral random matrix theories, from which the microscopic spectral correlation functions 
of the singular values can be determined. Our results can in principle be tested in future 
lattice QCD simulations. This should provide a value of the diquark condensate at any 
density, by which the conjectured BEC-BCS crossover could be confirmed numerically. 

It is evident from our results that the existence of a nonzero diquark condensate at 
any chemical potential implies an accumulation of the Dirac singular values at the origin at 
any quark density. In the case of the QCD vacuum, near-zero Dirac eigenvalues responsible 
for chiral symmetry breaking are believed to originate from instantons, as illustrated by 
the instanton liquid model [100]. One might thus naively expect that the accumulation 
of near-zero singular values in dense QCD is also attributable to instantons. However, 
this is not the case. Although instanton effects are presumably important at small and 
intermediate density, they are suppressed at sufficiently high density where the one-gluon 
exchange interaction is more important for the formation of diquark pairing. The presence 
of the Fermi surface is crucial in this mechanism. 

Let us discuss some possible extensions of the present work. First, here we have 
concentrated on nonzero diquark sources without quark masses, but the generalization to 
include quark masses looks straightforward. In that way we can study not only the diquark 
condensate but also the chiral condensate or the BCS gap as a function of the chemical 
potential. Second, the results obtained in sections 5 through 7 could also be generalized 
to theories with (3 = 2 and (3 = 4. Third, it would be very interesting to generalize 
our results to the color-superconducting phases of three-color QCD. Unfortunately, this 
is not straightforward since the diquark source is no longer gauge invariant, though the 
magnitude of the diquark condensate is. The object D(n)^ D(n) obviously exists in three- 
color QCD and can be studied, but it is currently unclear to us how its spectrum is related 
to physical observables. It would also be interesting to find out whether (and if so, how) 
the gauge invariant four-quark condensate is related to the Dirac eigenvalue or singular 
value spectrum in three-color QCD. 
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We conclude with a short discussion of a phenomenon analogous to the conjectured 
BEC-BCS crossover in QCD-like theories studied in this paper, i.e., the conjectured hadron- 
quark continuity in three-color QCD at nonzero baryon density [101]. Recently, the exis- 
tence of a bound state of baryons, the H dibaryon (with mass ms), was observed in lattice 
QCD simulations at zero density [102, 103]. This implies (for three degenerate flavors) that, 
as we go to nonzero baryon density, first BEC of H dibaryons occurs, since due to its bind- 
ing energy a bosonic H dibaryon has a smaller excitation energy m#/2 — hb (per baryon) 
than a baryon. In this BEC state U(l)s and chiral symmetry are broken spontaneously, 
and thus this state has the same symmetry-breaking pattern as the color-flavor-locked 
phase at high density, where both U(1)b and chiral symmetry are broken by the diquark 
condensate. Therefore these two phases can be continuously connected without any phase 
transition (hadron-quark continuity) as first conjectured by Schafer and Wilczek [101]. An 
explicit realization of this conjecture was given within generalized Ginzburg-Landau theory 
[104-106] and the Nambu-Jona-Lasinio model [39], where it was shown that the QCD axial 
anomaly can lead to a crossover between the hadronic and the CFL phase. Although it is 
not yet clear what happens physically between these two regions, there could be successive 
changes of states, from BEC of dibaryons to BCS pairing of dibaryons to BEC of diquarks 
and finally to BCS pairing of diquarks. 

In the real world, flavor symmetry is explicitly broken due to the heavy strange quark 
mass, and the existence of the H dibaryon as well as the scenario above may be modified. 
Actually, it is empirically believed that a nuclear liquid-gas phase transition takes place 
as the baryon density increases. In order to understand high-density matter in the real 
world, it would thus be crucial to take into account the effects of flavor symmetry breaking. 
Unfortunately, these effects cannot be directly studied in lattice QCD simulations even for 
QCD-like theories since, e.g., nondegenerate quark masses in two-color QCD give rise to 
the sign problem. It is an important future problem to figure out how one can study flavor 
symmetry breaking at nonzero density on the lattice. 
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A Definitions and conventions 

Unless stated otherwise we always work in Euclidean space. The 7-matrices are Hermitian 
and satisfy {7^, 7„} = 25^ (fi, v = 1, . . . , 4). We choose the chiral representation given by 

7 * ={ia i J' 74 = ^1 J ' * = ™™* = ^0-lJ' (A.l) 
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where the (Xj (i = 1, 2, 3) are the usual Pauli matrices. A Dirac spinor ?/; can be written in 
terms of two Weyl spinors i/jr and tpL, 

and we have 

^ = ^ 74 = tyt ^ (A _ 3) 
The projection operators on the right- and left-handed sectors are 

P K =i(l + 75) and P L = i(l- 75 ), (A.4) 
and we have, with a slight abuse of notation, 



(A.5) 



The charge conjugation matrix satisfies 

C 7/ ,Cr 1 = - 7 J and C T = -C. (A.6) 

We adopt the choice C = 27472, for which we have 

C- l = C ] = C, (A.7) 
[C, 71] = {C, 72} = [C, 73] = {C, 74} = [C, 75] = . (A.8) 

We also define the A'j-dimensional antisymmetric matrix 

I=( "M. (A.9) 

B Partition functions with diquark sources 

In this appendix we derive the singular value representations of the partition functions 
with diquark sources for the theories with /3 = 1, 2, and 4, taking exact zero modes into 
account and showing how the positivity of the path integral measure is determined. 

B.l Two-color QCD ((3 = 1) 

We first consider two-color QCD. We need to evaluate the Pfaffian of the operator W in 
(2.6) in the chiral limit, regarded as an infinite-dimensional antisymmetric matrix, 



Ct 2 (J r P r + J l P l ) —D([i) T 

D(jm) -Ct 2 (J r P l + JlP R ) 



Pf(W) = Pf (^v^-y^i,, „_, t w t „J. (B.l) 



Note that the transpose D(fi) T includes transposition of the space-time indices, in addition 
to the color and spinor indices. Below we will give two treatments of Pf (W). They lead to 
the same form for the partition function, which underscores the correctness of the result. 
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B.l.l Rigorous derivation 

To evaluate the Pfaffian it is useful to employ a specific functional basis. Here we use the 
eigenfunctions of D^D and DD^ introduced in section 3. 2, 48 

D^Dip n = £*(p n , J (fx iflnipn = 5 mn , (B.2a) 

DD ] (p n = £l<p n , J d A x filntfn = 5 mn , (B.2b) 

where for > we have 

ip n = ^D ] (p Tl and <p n = ^-D(p n . (B.3) 

Then the fields tp and ip in (2.4) can be expanded in the bases {y? n } and {^n}, respectively. 
We need the help of the following lemma to proceed. 

Lemma. Without loss of generality, we can assume 

= p n CT 2 ip* n , (B.4a) 

<Pn = -p n CT 2 ip* n (B.4b) 

for some p n G C with \p n \ = 1. 

Proof. For two-color QCD, it can be shown from (2.2a) that 

{D ] Df = Ct 2 D ] DCt 2 , (B.5a) 
(DD^)* = Ct 2 DD ] Ct 2 . (B.5b) 

Using these properties as well as (B.2), it follows that 

rtD{CT 2 y* n ) = e n {Cr 2 ^ n ) , (B.6a) 
DD^Cr 2 <p* n ) = e n (Cr 2 <f* n ). (B.6b) 

Therefore both tp n and CT 2 ip* n {(p n and Cr 2 ip^j are eigenfunctions of D^D (Z)Z)t) with the 
same eigenvalue Then two possibilities arise: 

1. ip n and Cr2(/9* are linearly independent, and the eigenvalue ^ is (at least) doubly 
degenerate. 

2. (p n and Cr 2 ip* n are linearly dependent, and the eigenvalue £^ is not degenerate. 

In the first case, we can redefine M\{(p n + p n CT 2 (p* n ) as ip n and N 2 (tp n — PnCT 2 (p* n ) as (p n +i 
with normalization constants A/12 an d arbitrary p n G U(l) so that 49 

= p n Cr 2 ip* n and ip n+ i = -p n Cr 2 ip* n+l = p n+1 CT 2 ip* n+1 , (B.7) 



48 Note that these eigenfunctions are ordinary c-number functions and not Grassmannian. 
49 This redefinition does not change the chirality because CV2 commutes with 75. 
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where we used (Ct 2 )* = Ct 2 and (CT2) 2 = 1. In the second case, it can easily be shown 
that there exists a phase p n 6 U(l) such that <p n = p n CT 2 tp* n , so the desired relation holds 
automatically. This completes the proof of (B.4a). Note that we have shown this for both 
zero and nonzero modes. 

Let us proceed to the proof of (B.4b). For nonzero modes (£ n > 0), we have 

= iz>X = U-Cr 2 DCr 2 )plCT2<Pn = -p* n Cr 2 ^Dcp n = -p* n Cr 2 ^ n . (B.8) 

Multiplying both sides by —p n Cr 2 we obtain (B.4b). Finally we consider the case of zero 
modes (£ n = 0). Although there is no obvious relation between the zero modes of D and 
those of D\ their numbers are equal, see (D.24). Choosing the phases of the zero modes 
of DD* one can make (B.4b) hold. □ 

Now it is straightforward to obtain the matrix elements of W in the bases {^p n } (for 
-0) and {(p\i} (for ijj). We consider separately the four blocks of W. 

• (l,l)-block: 

f d 4 X^l[CT 2 (J R P R + J L P L )]p n = f d A Xp m ^ m (J R P R + J L P L )Vn 

= Pn5mnJ R /L , (B.9) 

where in the first step we have used (B.4a) and in the second step R/L corresponds 
to the handedness of (p n . 

• (2,l)-block: 

J d 4 x<pl n Dip n = £ n 5 mn . (B.10) 

• (l,2)-block: 

J d 4 x <pl(-D T )v* n = -i n 5 mn . (B.ll) 

• (2,2)-block: 

J d 4 x (pi [ - Cr 2 (J R P L + jlP R )] <p* n = j d A x vl(J R P L + jlP R )Pn(pn 

= P*JmnJ R/L , (B.12) 

where in the first step we have used (B.4b) and in the second step R/L corresponds 
to the opposite of the handedness of (p n . 

Collecting all results, we obtain 

v ; 1 dmg(u) ® t Nf diag(K) ® 4 /L J A n A ^ e™ 4/ J v ; 
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where we used p n p n = 1- The contribution from zero modes can be read off as 

= [Pf(J R )Pf(j|)]" fl [Pf(4)Pf(J L )] ni , (B.14) 




where tir, l > denotes the number of zero modes of each handedness. 
The contribution from a nonzero mode is given by 50 

det l/2 I J R/L \ _ fdet 1/2 (^ + J R J R ) for right-handed <p n , 

\ & J R/L ) ~ \det 1 / 2 (^ + J{J L ) for left-handed <p n , !> " 

where we used the fact that the handedness of (p n is opposite to that of <p n for £ n ^ 0. Since 
(B.15) is manifestly positive definite, the (non-) positivity of the measure is determined by 
(B.14). 

Summarizing, we find that the full partition function of two-color QCD with diquark 
sources in the chiral limit reduces to the following expression in terms of the singular values 
ofD(p), 

Z{Jl,Jr) = /[Pf(JH)Pf(Ji)] na [Pf(4)Pf(J L )] Bi 



n'det 1 / 



/2 ' fn + JrJrPr + J&lPl ) ) , (B.16) 

YM 



where the primed product runs over all nonzero singular values. The above expression 
allows for accidental zero modes, but they are of measure zero, see section 3.1. Therefore 
generically we have (ur, til) = iy, 0) for v > and (0, —u) for v < 0. Note also that 

pf(4)pf(j L ) = [pf(j fl )pf(jt)]*. 

B.1.2 Short derivation 

The derivation of the singular value representation of the partition function in the last 
subsection is rigorous but lengthy. Here we will give a less rigorous but shorter derivation 
of the same expression. 

Let us assume /i = so that the extended flavor symmetry SU(2iVj) is intact in the 
absence of diquark sources. In computing Pf(W) let us separate the contributions from zero 
modes and those from nonzero modes. First, in the space of nonzero modes, we compute 
the square of the Pfaffian in (2.6), which is equal to its determinant, 

V D -Ct 2 (jIp l + j{p r ))' • ' 

where the prime on both sides indicates the omission of zero modes and we used —D T = 
Ct 2 D^Ct 2 in the (1,2) block. If we interchange the first and the second column, there will 
be a factor (— l) d , with d the total dimension of the space spanned by the nonzero modes 



50 Rigorously speaking, Pf and det 1 / 2 may differ by a sign. However, a ^-independent multiplicative 
constant can safely be omitted without changing expectation values of observables. 
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(we assume a suitable regularization such that d < oo). Since all nonzero modes are paired 
by 75, d is an even integer and (— l) d = 1. Thus 

Ct 2 D^Ct 2 Ct 2 (JrPr + JlPl)\ 
-Ct 2 {J r P l + jIp r ) d ) 

D^Ct 2 J r P r + J L P L \ 
-J R P L - J[P R Ct 2 D ) 





D^Cr 2 JrPr + J L P L 
rPl ~ 4i 



-AP L - J\Pr Ct 2 D 



(B.18) 



Using the formula det = det(AD — AC A l B) that holds when the blocks are 

square matrices of the same dimension and A is invertible, we obtain 

Pf'(W) 2 = det' (d*D + D^Ct 2 (J r P l + jIPr^Ct^-^JrPr + J L P L )) 

= det'(D^D + A R J R P R + A L J L P L ). (B.19) 

We stress that the above formula cannot be used in the full eigenspace of D because D^ 1 
does not exist if there are zero modes. Since (B.19) is positive definite we can take its 
square root naively to obtain a formula for Pf (W). 

The next task is to incorporate zero modes (where from now on we ignore accidental 
zero modes). We do this by switching from quark masses to diquark sources. Let us recall 
that the mass term in (2.3) can be written as [20] 

?(MP L + M^P R )iP = *#W 2 ^_ ° /T #* + h.c. (B.20) 

with the "extended" spinor Vl/ = ( R ] . 51 As is well known [4], the contribution of 

\a 2 T 2 ip L J 

zero modes to the fermion determinant in the absence of diquark sources is (det M^)" for 
v > and (det M)~ u for v < 0. 52 This expression is not desirable, as it does not manifestly 
show the spurious invariance under SU(2iVj). 53 Instead, we write detM as a Pfaffian [88, 
Eq. (4.13)], 



detM=(-iyV 2 Pf l T M \, (B.21) 




51 The symbol ^ used here is not to be confused with defined in (2.4). 

52 This differs from [4, 88]. The reason seems to be that our 75 differs from their 75 by a minus sign. (The 
definitions of the mass term and of v in (2.11) are identical.) Thus v — tir — til in this paper, whereas 
v = riL — riR in [4, 88]. One should be careful when comparing results in this paper with those in [4, 88]. 

53 For a fixed quark mass, SU(2iV/) is broken to U(l)s x SU(7V/)u x SU(AT/)i, but the initial symmetry 
is kept intact if we formally transform the quark mass as a spurion field, see (B.22). 
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which makes the invariance under U £ SU(2iVy) explicit, i.e., 






det U ■ Pf . ^ = Pf . ^ . (B.22 



Pf 

Similarly the diquark source in (2.3) can be cast in the form 

\VCt 2 (J r P r + J L P L )^ + h.c. = V<t 2 t 2 [~ 3 * °J + h.c. (B.23) 

Comparing this with (B.20) we notice that the diquark source is obtained if we simply 

replace ( ™ I by | 3r ^ | . Therefore the contribution of zero modes (for v < 0) 
y \-M T J \ J L J 

in the chiral limit and with nonzero diquark sources can be found by the replacement 54 
(detM)~ u = (-lf N f/ 2 Pf 




_^ { _ lT N f ,2 ¥i ( 4 0j = [Pf(jt)Pf(J L )]-", (B.24) 

and similarly for v > 0. Combining the square root of (B.19) with (B.24) we finally obtain 
for the partition function 

Z(J L , Jr) = ^ | [S f ( ™^ ( ^p } y/detfiDW + jynPR + jlJLPL)^ , (B.25) 

where the first (second) line in curly braces applies to v > (^ < 0). Assuming that 
/j / does not change the form of this expression, we arrive at (B.16) in the previous 
subsection. Note that this argument hinges on the assumption that the contribution of 
exact zero modes at fj, ^ is the same as at \i = 0. We also point out that, since D = —D* 
at n = 0, we could also have written D^D appearing in (B.25) as — D 2 , but this would not 
have generalized to /j / 0. 

B.2 QCD with isospin chemical potential (/3 = 2) 

Here we only sketch the shorter derivation of the singular value representation of the 
partition function, i.e., the extension of the argument in section B.1.2 to j3 = 2, and omit 
the lengthy and rigorous version. 

First, working along similar lines as in section B.1.2, the contribution from nonzero 
modes to the partition function can be shown to take the form 

det' {D^D + pp*P R + XX* Pl) . (B.26) 



54 A precise account of this procedure goes as follows. We start with a certain diquark source in the chiral 
limit. By a suitable SU(2JV/) transformation we can rotate it into the form of the mass term, for which 
we know that the zero modes contribute (detM) - " to the fermion determinant. Then we rotate inversely, 
bringing the mass term back to our original diquark source. This is how we get (B.24) so quickly. 
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Next we consider the zero modes. In order to figure out the mapping between quark 
masses and pionic sources, we go to a new basis 

••Hi-te) md *-4Hi)- <b - 27) 

These are the (3 = 2 counterparts of the extended spinor introduced after (B.20) for 
(3 = 1. Then the mass term for M = diag(m u , m<f) is 

u(m u P L + m* u P R )u + d(m d P L + m* d P R )d = * R MV L + V L M*V R . (B.28) 

The zero-mode contribution to the fermion determinant is given by (det M)~ v for v < and 
(detM^) 1 " for > 0. This form is manifestly invariant under the spurious SU(2)^ x SU(2)x 
rotation of M . On the other hand, the pionic source term reads 

u(X*P R + p*P L )d + d(pP R + XP L )u = V R * L + * L ^ X *^j V R . (B.29) 

Comparing with the mass term, we obtain the correspondence 

i.e., there is a mass term that can be rotated to a given pionic source term under the action 
of SU(2)^ x SU(2)i. Thus the contribution of zero modes is given by 

det (p 0* ) = (~ pX *^ fol "^>0, (B.31a) 

det ^ = (-p*\y u for i/ < . (B.31b) 
Summarizing, we find for the partition function 

A) = ^| det'C.Dt^ + ^p,, + AA^J^)^ for °J . (B.32) 

In deriving (B.32) we ignored accidental zero modes. The result can straightforwardly be 
extended to include them, and we then obtain (2.14). 

B.3 QCD with adjoint fermions (/3 = 4) 

Again we only sketch the shorter and less rigorous derivation. First, we consider the 
contribution of nonzero modes. Along similar lines as for (3 = 1 and 2 we obtain 



det'ptD + J R J R P R + J{J L P L ) , (B.33) 

but this is not the end of the story. Using (2.16a) we can show (D*D)* = CD^DC, from 
which it follows that a state Ctp^ associated with an eigenstate <p n of D*D is also an 
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eigenstate of D^D with the same eigenvalue. As ip n and Cc^* are linearly independent 
because of (CK) 2 = — 1, all eigenvalues of D^D are doubly degenerate (including zero 
modes). This allows us to take the square of (B.33), 

det"(£>t£> + JrJrPr + A J lPl) , (B.34) 

where in det" each degenerate eigenvalue of D^D is counted only once. 

Next we consider the zero- mode contributions. With the extended spinor ^ = 
of length 2Nf, 55 the mass term becomes V L 



MMP L + M*PM = ^ J *• - (_,° Mt -'fj * • (B.35) 

The contribution of the zero modes to the fermion determinant in the absence of diquark 
sources is given by [4, 88] 

det | ° * M | forI7<0 and det ( ° ~ %M | for V > , (B.36) 
\iM T J 1-iAft y ~ v 7 

where 1/ = (tir — n^)/2. 56 Because of the degeneracy mentioned above the numbers ur, 
of zero modes are even integers so that z7 £ Z. The number z7 is proportional to 
z/, the winding number of the gauge field. The proportionality constant depends on the 
color representation of the fermions. For example, V = N c v for fermions in the adjoint 
representation of SU(iV c ) [4]. More general cases are considered, e.g., in [107]. 
The diquark source term can be written as 

\VC( JrPr + JlPl)^ + h.c. = - 

(B.37) 

Comparison with (B.35) again suggests a correspondence between diquark matrix and mass 
matrix, from which the contribution of the zero modes in the chiral limit and with nonzero 
diquark sources is given by 








-J L 

Collecting everything, we find for the partition function 

rt 



det(- Jr J\) v for V > , (B.38a) 
det(- J r Jl)~ V for v < . (B.38b) 



det(-jU L ) 




Z(Jl, Jr) = ({ ^ J f J T L Zv \ det"(D^D + J R J R P R + j[j L P L ) 

' YM I 

B.39) 

Again, in deriving (B.39) we ignored accidental zero modes. The result can straightfor- 
wardly be extended to include them, and we then obtain (2.20). 

55 Again, the symbol ^ used here should not be confused with the ^ defined in (2.4) or after (B.20). 
56 Note that u in [4, 88] is — v in our notation, see footnote 52. 
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C Consequences of a (non-) positive measure 
C.l Diquark sources and positivity 

From the results of appendix B we can immediately read off conditions for the positivity of 
the measure. Here we consider some cases of particular physical interest (assuming 6 = 0). 

m /3 = 1: From (B.25) it follows that 

1. The choice Jr = —Jl (source for the + diquark condensate) implies 

pf(j fl )pf(4) = pf(j H )pf(-4) = (-i)^|pf(j fi )i 2 = |pf(ji?)i 2 > o. (c.i) 

Hence the measure is positive definite and there is no sign problem. 

2. The choice Jr = Jl (source for the 0~ diquark condensate) implies 

pf(j fl )pf(4) = pf(j fl )pf(4) = (-l^pfcjfl)! 2 . (c.2) 

This is positive (negative) if Nf = 4ra (Nf = 4n + 2) with n £ N. Thus the sign 
problem arises for Nf = An + 2 due to the topological sectors with odd v. 

• 13 = 2: From (B.32) it follows that 

1. The choice p = — A (source for the + pion condensate) implies 

-pX* = pp* > (C.3) 

so that the measure is positive definite and there is no sign problem. 

2. The choice p = A (source for the 0~ pion condensate) implies 

-p\* = -pp* < . (C.4) 

Thus the sign problem arises due to the topological sectors with odd v. 

• p = 4: From (B.39) it follows that 

1. The choice Jr = — Jl (source for the + diquark condensate) implies 

det(- J R J{) = det(J R J R ) > (C.5) 

so that the measure is positive definite and there is no sign problem. 

2. The choice Jr = Jl (source for the 0~ diquark condensate) implies 

det(- Jr4) = (-l) N f det(J*4) • (C-6) 

This is positive (negative) for even (odd) Nf, 57 so the sign problem arises for odd 
Nf due to the topological sectors with odd V. 

57 Note that Nf is the number of Dirac fermions (not Majorana fermions). 
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C.2 Sign problem and Dirac spectrum 
C.2.1 Discussion by Leutwyler and Smilga 

A long time ago it was pointed out by Leutwyler and Smilga [4, Sec. VII] that the positivity 
of the measure is a necessary condition for the Banks-Casher relation [1] to hold. They 
used Nf = 1 QCD at p, = as an example. Here we review their discussion briefly and 
motivate our microscopic analysis in the next subsection. 

A standard derivation of the Banks-Casher relation at /x = goes as follows. Assuming 
9 = for simplicity and taking m to be real (but not necessarily positive), we have 

kL i °s('" l " i n'«+'" 2 )> YM 

n 

1 „ ,, 1 /v^' 2m 

1 + TA (E 



m z / N f =i 

n 



where the i\ n are the Dirac eigenvalues, the primed sum stands for the summation over 
nonzero modes, and 

P(A) = ^(EV-A„)) (C.8) 

n ' 

Note that p(X) depends on m through the fermion determinant in the measure. The cases 
m > and m < differ qualitatively, as we shall see now. 

• m > 0: Since {V 2 )at /= i = mV^T, with S a low-energy constant, the first term in (C.7) 
is suppressed as 0(l/y/Vi) and becomes negligible at large volume. The second term 
in (C.7) will reduce to Trp(0) as m — > + after the thermodynamic limit. Therefore 

(^) Nf = 1 =7T P (0). (C.9) 

• m < 0: Since a probabilistic interpretation is no longer possible one cannot drop the 
first term of (C.7). Indeed we have [4] 

1 S e 2 l x l 
— — (\v\) Nf= i = -)■ oo as a; = ml^S ->• oo . (C.IO) 



/■°° 2m 

Similarly / d\ p(X) —5 =• diverges to —00, but these divergences cancel neatly 

in (C.7). 7o X +m 

Now it is clear why the Banks-Casher relation (C.9) fails for m < 0: 

• Zero- mode contributions are not suppressed at all in the thermodynamic limit. Non- 
trivial topologies cannot be ignored. 

• The macroscopic spectral density p(X) is ill-defined in the thermodynamic limit at 
least near the origin. 
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It seems natural to expect that these findings also apply to dense QCD-like theories, where 
the diquark (or pionic) sources can spoil the positivity of the measure (see section C.l). In 
a measure with indefinite sign, p sv (£) would be singular in the thermodynamic limit, and 
we cannot derive Banks-Casher-type relations for the diquark and pion condensate. This 
is the reason why in the main text we chose only those sources that respect the positivity 
of the measure. An exceptional case is the limit [i = oo, where instantons are completely 
suppressed and the exact zero modes will disappear. For this case we can in principle 
choose any diquark/pionic sources. 

We conclude this subsection with a few important remarks. 

• Sometimes there are multiple external fields that couple to different condensates and 
can be inserted without spoiling the positivity. For instance, in Nf = 4 QCD at \i = 0, 
the Banks-Casher relation can be derived for both (ififipf) and (V'/H'sV'/) because 
the degenerate purely imaginary mass term also respects positivity. 58 However, p(0) 
itself can be measured on the lattice without using any external sources. To which 
condensate p(0) is related to is determined by the infinitesimal external field we 
select. 59 

• The failure of a Banks-Casher relation does not prove the absence of the condensate 60 
because it is possible that the condensate exists but is not determined by the value 
of p(A) at the origin. An instructive example is three-color QCD at small [i. It 
supports a nonzero {iptp), but there is no Banks-Casher relation. The microscopic 
spectral density of the complex Dirac eigenvalues exhibits a drastic oscillation with 
a period ~ 0{1/V±) and an amplitude ~ exp(Vy [108, 109]. This suggests that the 
macroscopic spectral density in the limit V4 — > 00 is singular at least near the origin. 

In the discussion by Leutwyler and Smilga, a sick behavior of p(X) for Nf = 1 with 
m < was alleged in an indirect way. In the next subsection we study the spectral density 
in the microscopic limit (e-regime) and explicitly show how it behaves in an indefinite 
measure. 

C.2.2 The microscopic limit 

The effect of the topologically nontrivial sectors {y 7^ 0) on the Dirac spectrum was studied 
in detail in [110], where the microscopic spectral density and the chiral condensate in full 
QCD (including all topologies) were analyzed. All explicit examples in that reference 
were worked out for positive quark mass and vanishing CP-breaking angle 8 = 0. Here we 
demonstrate the effect of negative mass, or equivalently nonzero 9, on the microscopic Dirac 
spectrum. As in [110], we use a formula [111] that expresses the microscopic correlation 
functions in terms of finite-volume partition functions. While in [110] the contributions 
from different topological sectors were simply summed up numerically, we here use a closed 

58 For Nf = 2 this is not the case. This fact does not seem to have been appreciated in [50]. 
59 More precisely, we have | (ipip) | = | (^175^) | = np(0), and the orientation of the condensate is deter- 
mined by the external field. See also footnote 24. 
60 Indeed (qq)Nf=i is the same for m > and m < 0. 
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Figure 8. Microscopic spectral density of the Dirac eigenvalues in N c > 3 QCD (j3 = 2) for 
Nf = 1 at 9 = 7r for different quark masses. Here p s (x) is not necessarily positive definite and the 
probabilistic meaning as a "density" of eigenvalues is lost, reflecting the underlying sign problem. 

expression for the full partition function summed over all topologies that became available 
later [112]. By doing this we avoid the errors coming from the truncation of an infinite to 
a finite sum. 

Since all relevant formulas can be found in the literature [110-112] we omit the technical 
details and only show the final plot. As the simplest case we considered Nf = 1 QCD at 
p = with N c > 3 (p = 2). We have chosen 8 = ir and m > so that the measure is not 
positive definite. (This is physically equivalent to 9 = with m < 0.) In figure 8 we show 
the plot of the microscopic spectral density with several values of the quark mass. In the 
chiral limit (m = 0) the contribution from the v = sector is dominant and there is no effect 
of 8. For nonzero quark mass, the microscopic spectral density exhibits a strong oscillation 
whose magnitude grows rapidly with mV^E, while the period of oscillation looks roughly 
independent of mV^T,. This is exactly the opposite of what we observe at 8 = 0, where 
the dynamical quark decouples for a sufficiently large mass and p s (x) becomes smoother. 
This oscillatory behavior and the failure of decoupling is reminiscent of the complex Dirac 
eigenvalue spectrum in QCD at p ^ [108]. We also checked that these characteristics are 
not limited to 8 = tt: Even for a small 9^0 the violent oscillation shows up eventually at 
a sufficiently large value of mV^T,. 

Now we have understood microscopically why the standard Banks-Casher relation can 
fail: The spectral density varies rapidly over the scale of the quark mass and p(0) is 
certainly ill-defined in the thermodynamic limit. It is quite tempting to conjecture that 
this phenomenon is universal and will occur also in other theories with indefinite measure. 
In particular, we expect this for dense QCD-like theories with positivity-breaking external 
sources. 61 Indeed, in Nf = 2 two-color QCD, jn = ji at 9 = is equivalent to jr = —jz 
at 8 = 7T, and the similarity to the case considered above is evident. It would be interesting 

61 The ^-dependence of QCD-like theories at fi ^ was studied in [113, 114], but the singular value 
spectrum of the Dirac operator was not considered there. 
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to study the impact of the sign problem on the singular value spectrum in more detail in 
analogy to the analysis of [108, 109]. 

C.3 QCD inequalities 

QCD inequalities are a powerful theoretical tool to impose strong constraints on the dy- 
namics when the Euclidean functional measure is positive definite [115-118]. Although the 
sign problem hampers the application of QCD inequalities to three-color QCD at // 7^ 0, 
they are still applicable to a class of dense Q CD-like theories with positive definite measures 
[19, 21, 22, 26]. In particular it was shown for two-color QCD (with nonzero quark masses 
and no diquark sources) that the symmetry breaking is driven by the + diquark conden- 
sate [19], thus leaving parity unbroken. In this appendix we reexamine the application of 
QCD inequalities to two-color QCD at /j, 7^ and extend the analysis to the case of nonzero 
diquark sources. For the sake of simplicity we will limit ourselves to either nonzero quark 
masses or nonzero diquark sources, but not both. We always work in Euclidean space. 

C.3.1 Two-color QCD with quark masses 

We review the original discussion in [19] and assume even Nf for positivity of the measure. 
In this subsection, the transpose (T) and the adjoint (f ) act only on color, spinor, and flavor 
indices, but not on space-time indices. The quark chemical potential jj, can be arbitrary 
in the following. When the diquark sources are absent and the (degenerate) quark mass is 
real and positive, the propagator in a fixed background gauge field is given by 

where ( )w, denotes the average only over the fermion fields. Below we also use ( )^ t A an d 
( )a for the full average and the average only over the gauge fields, respectively. 

Let us consider a diquark operator M(x) = ip T Fip with an antisymmetric matrix F 
satisfying T^T = 1. The correlation function of this field is then given by 

(M(x)M*(y)) = U T {x)T^{x)^(y)Tij T {y)) (f = 74^74) 

= (tr [TS^(x,y)TS^(x,y) T ]) A 

<(tr[S^(x,y)S^(x,y^]) A , (C.12) 

where the last line follows from the Cauchy-Schwarz inequality and the positivity of the 
measure. The same upper bound can be proven for those correlators of mesonic fields ipTtp 
that have no disconnected piece. 62 On the other hand, (D + m) T = Ct2^(D^ + m)Cr275 

62 By "disconnected piece" we mean a diagram whose quark lines are not connected directly, although 
they could be connected by gluons, see, e.g., figure 1 in [116]. The absence of disconnected pieces is essential 
for QCD inequalities. The correlators of a and rf contain disconnected pieces and our arguments do not 
apply. Indeed, if the disconnected piece could be dropped, the rj' would be as light as the pions. This 
is of course invalid. In reality the contribution of the gluonic intermediate state \FF) is essential in the 
generation of the large mass of the 77'. 
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allows us to write 

(M(x)Mt(y))^ = (tr [TS^(x, y)TCr 2l5 S^(x, y)^Cr 2lh ] ) ^ . (C.13) 
If we assume V = Ct2^^T j for an antisymmetric flavor matrix Tf, it follows that 

(M(x)Mt(y))^ = tr[r>rt] (tr [S^x,y)S^{x,y?] )^ . (C.14) 

This is equal to the upper bound in (C.12) up to an irrelevant multiplicative constant, 63 
implying that the diquarks in this channel are the lightest of all mesons and diquarks. 64 
For T = CT2T f the upper bound is not saturated, implying that if a diquark condensate 
forms, it does so in the + channel and not in the 0~ channel. This is consistent with 
the observation that the instanton-induced interaction for Nf = 2 is attractive in the + 
diquark channel and repulsive in the 0~ diquark channel [69, 119, 120]. 

It is a general rule that the preferred direction of a condensate depends on the direction 
of the external symmetry-breaking field. In the case above, we considered real positive 
quark masses, but it is also instructive to consider other cases. Let us take degenerate 
purely imaginary masses for instance. For the sake of positivity (det M > 0) we require 
Nf = An with n 6 N. The propagator now reads 

and because of {D + im^) T = —Ct2(D + 277175) 'Ct2 the propagator satisfies S,-^(x,y) = 
— Ct2S^z(x, yyCr2- Thus we find from the second line of (C.12) that 

(M{x)M\y))^ A = - (tr [TS^(x, y)TCr 2 S^(x, y^Cv,] )^ . (C.16) 

When T = Ct2Tj, the r.h.s. reduces to tr[rjTj]^ tr [S^(x, y)S^(x, y)t] saturating the 
inequality (C.12). Therefore this time the diquark condensation occurs in the 0~ channel, 
breaking parity. The conclusion is that the quantum numbers of the condensate are not 
entirely determined by the internal dynamics of the system, but are sensitive to the external 
symmetry-breaking fields. 

The reader may worry that this is in conflict with the Vafa-Witten theorem [121] stating 
that parity is not spontaneously broken in vector- like theories. Let us make two remarks 
on this point. First, various loopholes in the original "proof" have been discussed in the 
literature [122-128], and therefore it cannot be regarded as an established mathematical 
theorem. Second, our result does not contradict recent work [128] on a Vafa-Witten-type 
theorem for fermion bilinears, since the positivity of the probability distribution function 

63 The mass is given by the exponential decay of the correlation function, for which the multiplicative 
constant is irrelevant. 

64 Strictly speaking, this argument does not exclude the possibility that there are other mesons or diquarks 
of the same mass. For example, at fi = there are mesonic states that have the same mass as the scalar 
diquarks due to the extended flavor symmetry. This degeneracy is lifted at \jl 7^ 0. 
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of the observable, which plays an essential role in [128], is not ensured for purely imaginary 
masses. 

We end with a brief remark on the rotation of the condensate. It is known from chiral 
perturbation theory that the chiral condensate at fi = rotates into the diquark condensate 
for fj, > m n /2. For real positive mass, {ipip) at [x = rotates into the + diquark condensate 
(ifjC'ysifj) [20]. For imaginary mass (and Nf = An), (ifti'ysip) at fj, = rotates into the 0~ 
diquark condensate (ipCip).® 5 

C.3.2 Two-color QCD with diquark sources 

We now turn to two-color QCD in the chiral limit with diquark sources. In this subsection, 
the transpose and the adjoint again act only on color, spinor, and flavor indices, with one 
exception: as in the other parts of this paper the adjoint in acts on all indices. 

In calculations of correlators, it is important to know whether there is a disconnected 
piece or not. If there is, the contribution of the gluonic intermediate states will invalidate 
the derivation of QCD inequalities. In the previous subsection, the annihilation of the 
diquark into gluons was trivially prohibited by the \J(1)b charge conservation, but the latter 
is explicitly violated once we insert diquark sources. Indeed, in the cases of Jr = ±Jl = jl 
considered below, the correlators of ^ T Cr2l'^ and ijj T Ct2^Ii^ have disconnected pieces, 
and QCD inequalities do not apply. We must limit ourselves to those correlators that have 
no disconnected piece. 

The propagators in a fixed gauge field are given by the inverse of W defined in (2.6), 

(S^y) y*.y)\ s (ffi™* <™M _ (c , 7) 

The operator inverse W -1 depends on the choice of the diquark sources. We first discuss 
the diquark sources with positive parity and then those with negative parity. As stated 
above we work in the chiral limit. 

For the + diquark source, i.e., Jr = — Jl = jl with real j, we find 66 

1 „ 1 

(C.18) 



The positivity of the measure is ensured for even Nf (see section C.l). We now consider 
the correlator of the diquark field M(x) = t/) T rt/), assuming that there is no disconnected 
piece. A rerun of the arguments leading to (C.12) then shows that 

M{x)M\y)) = U T (x)T^(x)^(y)T^ T (y)) (F = 74^74) 



65 More general cases, not necessarily preserving the positivity of the measure, were addressed in [113] 
based on chiral perturbation theory with a nonzero #-angle. 

66 Note that the denominator in (W~ 1 )22 contains DD t and not D' D. Using W in (2.6) one can explicitly 
confirm WW' 1 = 1. 
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< ( tr [S^{.x,y)S^(x,y)^ A . 
We now use (C.18) to show that 67 



(C.19) 
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D 



j 2 + DW 
1 



p + DW 
= CT 2 S^(y,x)CT 2 . 

Using this relation in the second line of (C.19) we obtain 



M{x)M\y) 



i>,A 



It is easily shown that for both T = CT 2 Tf and T = CT 2 ^5Tf, we find 



M{x)M\y))^ A = tr[r>rt] (tr [S^(x, y)S^(x, y)t] ^ 



(C.20) 



(C.21) 



(C.22) 



which saturates the inequality (C.19) up to a multiplicative constant. Thus i(j t Ct2T ftp 
and tp T Ct2^V ftp are the lightest diquarks (with degenerate masses). 

From oc I we see that the nonexistence of disconnected pieces is ensured if 
tr[iTj] = 0. Recalling from section 6.1 that an Nf x Nf antisymmetric matrix can be ex- 
panded in the generators t A of XJ(Nf) /Sp(Nf), we find that this condition is satisfied for all 
but not for A = 0. Thus we conclude that {tp T CT2t A tp}A^o and {ip T CT2j5t A ip}A^o 
are the lightest diquarks (with degenerate masses), whereas no information on tp T Cr2ltp 
and ip T Ct2^I'4> can be gained from QCD inequalities. 

The fields {tp T CT2t A tp} a^o an d {V ;T C' r 275^ 4 V ; }A^o can be interpreted as the NG modes 
associated with the spontaneous symmetry breaking SU(iVy)# x SV{Nf)i — > Sp(iVy)# x 
Sp(Nf)L. However, this breaking is caused by both {tp T Cr2ltp) and {tp T CT2')bIip) , which 
means that our result does not yield any constraint on the parity of the condensate. (This 
is to be contrasted with QCD with quark masses and no diquark sources, where we could 
determine the parity of the condensate even though we could not apply QCD inequalities 
to the correlators of a and 77'.) 

We now turn to the 0~ diquark source. For Jr = Jl = jl with real j we find 
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Ct 2 D 



DW 
1 



jCr 2 I 



j 2 + DW 



j 2 + DW 



Cr 2 -jCr 2 I 



(C.23) 



j 2 + DD\ 



For the positivity of the measure we require Nf = 4n (see section C.l). Since the off- 
diagonal entries are not changed from (C.18), our previous discussion goes through with- 
out modifications, and therefore we reach the same conclusion that {tp T CT2t A tp}A^o an d 



67 Note that the interchange of x in y in the second line of (C.20) is correct even though the outer adjoint 
does not act on space-time indices. 
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{V> T Ct275^V'} j 4^o are the lightest diquarks. Again, the parity of the condensate cannot 
be determined from this argument alone. 

A further discussion of the parity of the diquark condensate is given in the last para- 
graph of section 5.1. 

D Anomaly and index theorem for non-Hermitian Dirac operators 

In this appendix we prove the following properties of the Euclidean Dirac operator 68 




D(ji) = 7 v D v + ijrr4= I " (D.l) 



at zero temperature and nonzero quark chemical potential fi in a given gauge field: 

1. The chiral anomaly equation in the chiral limit, 

iNf 

with J^ u = itpj u j5ip, is unchanged. 

2. We have 

dim ker Dr — dim ker D^ R = v , (D.3a) 
dim ker — dim ker D^ L = —u , (D.3b) 

where v is the winding number of the gauge field, defined in (2.11). 

3. For /x = 0, (D.3) reduces to the ordinary index theorem 69 

dimkerZ)/? — dim ker Dl = v (D.4) 

since Dr = —D' L for \x = 0. 

The in-medium chiral anomaly and its impact on the phenomenology have been dis- 
cussed in numerous studies [39, 54-57, 69, 104, 119, 120, 129-140]. So far the obser- 
vation (D.2) has been made via perturbative calculations of the triangle diagram [130— 
132, 135, 138] as well as by Fujikawa's path integral method [131, 135, 138]. The latter 
encounters some subtleties at [i ^ that were not stressed in the preceding works. Also, 
in contrast to the chiral anomaly equation, the modification of the index theorem due to 
/i ^ was rarely considered so far. 

In this appendix we take a closer look at these issues. After some preliminaries in 
section D.l, we present in section D.2 two derivations of the anomaly equation (D.2) via 

68 The color generators are assumed to be in the fundamental representation of SU(iV c ) for N c > 2. 

69 Whether the r.h.s. of this equation is v or —v is determined by the convention for 75. Here we use 
75 = 71727374 (see appendix A), which leads to +u. Some authors use the convention 75 = —71727374, 
which leads to —v. As mentioned in footnote 52, it seems that in [4] (although not stated explicitly) the 
second convention is used, which differs from ours by a sign. 
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the path integral method in a style that differs from [131, 135, 138]. In section D.3 we then 
study the index of the Dirac operator at fi ^ 0. We will show that (D.3) can be proven for 
generic gauge fields and fx ^ 0, while (D.4) requires an additional input. In section D.3. 3 
we will discuss where the difference between (D.3) and (D.4) comes from. 

Three caveats are in order. First, we note that Fujikawa's analysis in its original form is 
equivalent to one-loop perturbation theory [141-144] although it is sometimes incorrectly 
said that it offers a nonperturbative derivation of the anomaly. Our analysis should be 
seen as the extension of this equivalence to \i ^ 0. Second, our analysis does not extend to 
lattice fermions directly. It is beyond the scope of this paper to analyze the effects of fi in 
the lattice formulation (see e.g., [145-147]). Third, the entire discussion in this appendix 
will be given for a fixed gauge field background. 

D.l Preliminaries 

D.l.l Path integral measure 

To clarify our stance on the subject and fix the notation, we begin by reviewing Fujikawa's 
method at \x = [148, 149]. It starts with the expansion of the fermion fields ip and ip in 
terms of the eigenstates {V'n} of D(0), 

= a n ip n (x) = ^2 a n{x\n) , (D.5a) 

n n 

^(x) = ^6„Vt(x) = ^6 n (n|x), (D.5b) 

n n 

where the a n and b n are Grassmann variables. The transformations from tp and ip to {a n } 
and {bn} are unitary by virtue of the orthonormality and completeness of {ip n }, which are 
guaranteed since -D(O) is anti-Hermitian. The functional measure then transforms as 

\\Vil)(x)V4>(x) = (det(n|x)) _1 (det(x|m))" 1 J|da n d6 n 

x n 

= (det(n|m)) _1 J^J da n db n = da n db n , (D.6) 

n n 

where det(n|x) and det(x|m) each represent the determinant of an infinite-dimensional 
matrix specified by the indices n, m, and x. By standard calculations the evaluation of 
the Jacobian for an infinitesimal chiral transformation leads to an infinite sum, tr7s = 
Yin ^nix)^ l54>n(x) , which is to be regularized by, e.g., e° 2 ^ 2 with A — > oo. 

Two remarks are in order. First, for the change of variables (D.6) to be unitary, the 
operator whose eigenstates are used to expand the fermion fields must be (anti-) Hermitian. 
Second, the final result does depend on the choice of the operator. Had we used the plane- 
wave basis, the result would simply have vanished: tr(7s e® 2 /^ 2 ) = 0. This is sometimes 
ascribed to the lack of gauge invariance, but the gauge invariance is merely a necessary 
condition. Indeed there is an instructive example, tr(7s q d ^ d "/^ 2 ^ = 0, which reveals that 
a seemingly unitary transformation from the gauge invariant eigenspace of D U D U to that 
of D is actually non-unitary in the presence of regularization [150]. 
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The bottom line is that the operator used to define the functional space must be 
chosen so as to diagonalize the fermion action [78, 149]. The operators D U D U and d 2 are 
not eligible, because their eigenbases do not diagonalize the action. In the following, we will 
work with this criterion as a guiding principle. This completes our preliminary comments 
on the chiral anomaly at fi = 0. 

D.1.2 Remarks on the literature 

D{(j) is no longer anti-Hermitian at \x ^ and the eigenstates lose orthogonality, spoiling 
the unitary transformation (D.6). This requires us to choose a more appropriate definition 
of the path integral measure. 

• In [131] the fermion fields were expanded in eigenstates of D(ifx), the Dirac operator 
with imaginary chemical potential. With D(ijj) being Hermitian (in their conven- 
tion), the eigenstates are orthonormal and the anomaly equation follows straightfor- 
wardly. 

• In [135] the fermion fields were expanded in eigenstates of D(0). Again the derivation 
of the anomaly equation is straightforward. 

• In [138] the fermion fields were expanded in right and left eigenstates of D(fi), but the 
subtlety is that the "eigenstates" were simply defined by multiplying the eigenstates 
of -D(O) by e ±/iX4 . They either blow up exponentially at infinity for T = 0, or spoil the 
anti-periodic boundary condition in the 24 direction for T > 0. Note also that this 
multiplication does not change the eigenvalues from purely imaginary ones at \x = 0. 
We suspect that the eigenstates thus obtained are not the legitimate eigenstates of 
D(fi) (see, e.g., [56, 57]) since the contribution from occupied states below the Fermi 
surface is not taken into account [57]. 

The schemes based on D(i/j,) and -D(O) [131, 135] do not meet our criterion, as the eigen- 
states of D(ifj>) and -D(O) do not diagonalize the action. While the final result for the 
anomaly equation obtained in [131, 135, 138] is correct, the rigorousness of these approaches 
is not completely obvious to us. 

D.2 Proofs of (D.2): Path integral derivation of the anomaly at /1 7^ 

In this subsection we will present new derivations of the anomaly equation (D.2) by the path 
integral method, which we believe place the results in the literature on a firmer footing. 
For simplicity we assume Nf = 1. The extension to Nf > 1 should be straightforward. 

D.2.1 Derivation based on D^D 

A non-Hermitian Dirac operator occurs not only in QCD at jx 7^ but also in chiral gauge 
theories where the gauge interaction involves a 75-coupling. We shall apply one of the 
methods devised for this situation [151] 71 to QCD at \i 7^ 0. We use the same orthonormal 

70 Note also that the index theorem derived in [138] corresponds to (D.4), while the correct result at fi 7^ 
is (D.3) if there are accidental zero modes. 

71 In chiral gauge theories this regularization is known to give the so-called covariant anomaly [148]. 
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bases {<p n } and {(f>n} as in (B.2). Note again that the £ n are not the eigenvalues but rather 
the singular values of D, with £ n > for all n. The orthogonality and completeness of 
these bases follow from the Hermiticity of D^D and DD' . The fields can be expanded as 

ip(x) = ^2a n ip n (x) , ip(x) = ^bnfilix) , (D.7) 

n n 

and the fermionic part of the action is then diagonalized, 

S = f (fx ^(x)D^(x) = ZnKa n ■ (D.8) 

Therefore this scheme meets the criterion mentioned above. Note that, had we expanded 
both i/j and ip in only one of the bases, the action would not have been diagonalized. 

Under an infinitesimal chiral transformation ip — > e ia,y5 ip, ip — > ip e ia75 , the action 
changes as 

e~ s e~ s ' = exp - S + J (fx a(x)d l/ J 5u (x) , (D.9) 
while the functional measure also changes as 

VipVip -)■ VipVip expj-i / d 4 xa(x) ^ ^(x)75V? n (x) + ^(x)75<^ ri ,(x)J | 

= lim VipVip exp | — i d 4 x a(x)^2 ^Pn( x )l5 e~ D ^ D l^ <p n {x) 

+ ^t ( ^)75e- DDt / A2 ^(x)]}, (D.IO) 

where a Gaussian cutoff has been employed. Because the phase factor a(x) is arbitrary, 
we have from (D.9) and (D.IO) 

d u J 5u = i lim V U (x) 75 e~ DtD / A2 <p n {x) + <^ (*) 7 5 e~ DD '/ A2 (p n { x ) 

n 

where the last line is obtained using standard algebra involving the plane-wave basis. The 
//-dependent terms appear in the calculation but disappear from the final result. We have 
thus obtained the conventional anomaly equation (D.2). 

D.2.2 Derivation based on D 

In the second derivation, we apply the method of [152, 153] to QCD at fi / 0. 72 Let us 
denote the right and left eigenvectors of D(fi) by {ip n } and {x-n}, 

D(n)ip n = Xn^n , xLDfa) = X n xi • (D.12) 



72 In chiral gauge theories this regularization is known to give the so-called consistent anomaly satisfying 
the Wess-Zumino consistency condition [148]. 
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The eigenvalues A n are no longer purely imaginary but complex in general. We assume 73 
that they satisfy the following orthonormality and completeness relations, 74 

/ d 4 X xln(x)lp n {x) = 5mn , ^ lpn(x)xl(y) = S*(x - y) . (D.13) 

J n 

The orthonormalization is possible because the left and right eigenstates corresponding to 
different eigenvalues are orthogonal, as can be seen from 



\ m J d 4 XxLi>n = J d A x{Xm D )^n = J d A Xx\n{ D ^n) = An / d 4 X xln^n ■ (D.14) 

Let us expand the fields ifi and ifi in the sets {ip n } and {xl}, respectively, 

H x ) = ^CLn^nix), lf)(x) = ^^nXn^) , (D.15) 

n n 

by which the fermion action is diagonalized as desired, 

S = J d 4 xipDip = ^2^2 J d A xb k x\{x)Daiil)t>{x) = ^ A fe 6 fc a fc . (D.16) 



k £ 



Thus this choice of bases also satisfies our criterion. After some algebra, the change of the 
measure under an infinitesimal chiral rotation is found to be 



VtpVip -> VipV^p exp 



lim P^P^exp 

A— >oo 



d 4 x 2a(x) ^2 xi(x)j 5 ^n(^ 

n 

■i J d A x 2a(x) ^ xi(x)l5 1 _ ^ 2/ / A 2 • ( D - 17 ) 



We avoided a Gaussian cutoff because complex eigenvalues are not suppressed by a Gaus- 
sian factor. Using the assumed completeness (D.13) we can move to the plane- wave basis. 
After some standard algebra, we find that the /i-dependent terms disappear from the final 
result and recover the anomaly equation (D.2). 

D.3 Proofs of (D.3): Index theorem at fi ^ 

In the following we present two derivations of the index theorem at \i ^ 0. The first 
derivation yields (D.3), while the second one yields (D.4). We analyze the origin of this 
discrepancy and relate it to the (in-) completeness of the bases. 

D.3.1 Derivation based on D^D 

Let us reexamine the analysis in section D.2.1. Prom (D.ll) it follows that 

lim V f d A x U(x) 75 e~ D '°/ A2 v n (x) + 4(x) 7 5 q~dd^ /A 2 -^1 = 2v (D lg) 



73 This assumption can be violated on a gauge field set of measure zero, see section D.3. 3, and then the 
derivation no longer works. 

74 The conditions in (D.13) say nothing about inner products such as J d 4 x ipj n (x)tp n (x). 
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Using (B.3) one can show that the entire contribution to the l.h.s. comes solely from the 
eigenstates with £ n = 0, for any finite A. This observation, combined with 

J d 4 x^ n (x)j 5 e- DtD/A2 ^ n (x) = dimkerL> i j-dimkerD L (D.19) 



n:£„=0 



and 



J d 4 x^ n (x)j 5 e- DDt/A2 ^ n (x) = dimkerL>[-dimker J £>J ? , (D.20) 

n: £n=0 

proves the equality 

dimker Dr — dimker Dl + dimker D^ L — dimker D R = 2v . (D-21) 

On the other hand, using the identity 

= tr(e- DZ)t /A 2 _ e -^/A^ (D22) 

that follows from the cyclic invariance of the trace 75 and the fact that all nonzero eigen- 
values of D^D and DD^ coincide, we have 

dim ker D = dim ker , (D.23) 

or equivalently 

dim ker Dr + dim ker Dl = dim ker D R + dim ker D^ L . (D.24) 

Then (D.21) and (D.24) prove (D.3). 

For v > 0, dimkerDft > v follows from (D.3a). For v < 0, dim ker Dl > \v\ = —v 
follows from (D.3b). Thus our index theorem at [i ^ shows that D(fi) must possess at 
least zero modes for any /i. This is consistent with the existence of exact zero modes 
in the instanton background at [i ^ [55, 56]. 76 

D.3. 2 Derivation based on D 

Let us return to the second derivation of the anomaly in section D.2.2. From (D.17) and 
the anomaly equation (D.2) it follows that 

If A n ^ 0, we have 

75 It was pointed out in [154] that the cyclic invariance of the trace could potentially break down. This 
does not happen in the instanton background, see [154] or section 4.2 of [155]. This possibility does not 
invalidate our result (D.3) since we have also constructed an alternative proof of (D.3), using the eigenbases 
of D${n) and Ds(— /i), which does not make use of (D.22). For brevity we do not show this proof here. 

76 Note that the fermionic zero modes in the instanton background at fj, ^ have a norm that diverges 
logarithmically in the spatial volume [56]. 
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= ~X£ 0*075 1 _ X 2/ A 2 D ^( X ) 

= -AnXn(^)75 1 _ ^rn^ ^nix) , (D.26) 

and hence all contributions to the l.h.s. of (D.25) from nonzero modes vanish. Restricting 
the sum in (D.25) to zero modes, we arrive at 



? 



dim ker Dr — dim ker Dl = v . (D.27) 

This "proves" (D.4), i.e., the index theorem in the same form as at \i = 0. It does not 
imply (D.3) proved in section D.3.1. The origin of this discrepancy is analyzed below. 

D.3.3 Analysis of the discrepancy 

The arguments in section D.3.1 proved (D.3), whereas in section D.3. 2 we were led to (D.4). 
Which is the true index theorem at fj, 7^ 0? 

To answer this question, let us begin by recalling a standard argument at \x = about 
the stability of the index of the Dirac operator against small perturbations. Suppose there 
are v right-handed zero modes and no left-handed zero modes in a given background gauge 
field. Now we deform the gauge field smoothly so that two of the eigenvalues A and —A 
(with eigenstates ip and 75^, respectively) approach zero. When A = is achieved, we 
have two accidental zero modes, which we can arrange to be eigenstates of 75 by choosing 
the linear combinations (1 + 75)^ and (1 — 75)V ; - Thus the number of right-handed (left- 
handed) zero modes becomes v + 1 (1), but their difference never changes: (y + 1) — 1 = v. 
This is why accidental zero modes cannot change the index. 

This argument is correct at \x = 0. If it could be extended to fx 7^ 0, the accidental zero 
modes could not change indD = dim ker Dr — dim ker Dl and the ordinary index theorem 
(D.4) would continue to hold. However, the argument does not work for fi ^ 0. The pitfall 
is that, when D(/j,) is not anti-Hermitian, ip and 75^ may become linearly dependent in 
the limit A — ¥ 0. Let us exemplify this by a toy model which mimics D(fi): 



Aoy(a) 



/0 2 -4\ 
12 
12 

\0a0 ) 



with eigenvalues {±2, ±V2a} . (D.28) 



The right eigenvectors of D^ y{a) associated with the eigenvalues v2a and — y2a are 
(4, — 2, 0, — \^2a) T and (4, — 2, 0, \^2a) T , respectively. In the limit a — > they coincide, 
implying that dim ker Dr increases by 1 whereas dim ker Dl remains zero. The accidental 
zero modes changed the index. 

At a = the dimension of the eigenspace of Dt oy (a) is 3 (< 4). This means that the 
eigenvectors fail to form a complete basis, a phenomenon that occurs only for non-Hermitian 
matrices. The interpretation of the discrepancy between the last two subsections is now 
straightforward: Our derivation in section D.3. 2 led to (D.4) because we assumed the 
completeness of the basis in (D.13). If the completeness is not ensured, it is not possible 
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to expand an arbitrary fermion field in this basis as in (D.7), and then our discussions in 
section D.2.2 and D.3.2 break down. 77 The index theorem that holds for generic gauge 
fields and /i / is not (D.4), but (D.3). 

However, the incompleteness of the basis requires a fine-tuning of the gauge field, and 
so we expect that this occurs only on a gauge field set of measure zero. For practical 
calculations it is justified to neglect the possibility of incompleteness, in the same sense as 
we can neglect accidental (non-topological) zero modes at \x = 0. Thus both (D.3) and 
(D.4) will be valid for almost all gauge fields and /i 7^ 0. Summarizing, we have 

for v > : dimker Dr a = dimkerl?^ a = v , dimker D R =' dimkerZ)^ =' , (D.29a) 
for v < : dimker D R a = dimkerD^ a =' — v , dimker Dr a =' dimker D^ L =' , (D.29b) 

where =' denotes an equality that holds "almost surely" . This completes our discussion of 
the index theorem at fi 7^ 0. 



E Derivation of (6.19) 

In this appendix we outline the derivation of the singular value density in the free limit in 
(6.19). The basic relation we use is 



where e — > + is tacitly assumed. To obtain the resolvent on the r.h.s., we express it in 
terms of the free Dirac operator D(ji) = ( r y v d v + ^a)®^-n c (/■* > 0) and move to momentum 
space, 

&e n -e+ie) = ( tr wmfi -e+ie) (K2) 



d A v 1 

V*N C I — -Vrtr- - -= . (E.3) 

' (2^) 4 {-ii> + ^i){if + ^i)-e + ie 



Then 



1- (£ = ~^ v ^ Jj2^\p 2+ » 2 -?\ 6 W + » 2 ~ ^ 2 ~ 4 ^ 2 p 2 ) • 

(E.4) 

where p = (pi,P2,P3) denotes the three-momentum. After elementary but tedious algebra 
this integral can be performed analytically and (6.19) follows. 



77 Although our discussion here is based on finite-dimensional matrices, we believe that the essential part 
of the argument carries over to the actual Dirac operator as an infinite-dimensional matrix. 
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F Random matrix theory for QCD with isospin chemical potential 
(0 = 2) 



In this appendix, we comment on the random matrix theory describing QCD with three 
colors, two flavors, and isospin chemical potential pi = 2p, see sections 2.2 and 4.2. Starting 
from the two-matrix model of [156] the random matrix theory for this case in the chiral 
limit is 

Z™ Nf = 2 ^h,h) = J dAdB e-^AU+Bt^det ( D W * J (F.l) 



with 

where ^4 and B are complex iV x (iV + v) matrices and 

X=(^ N ° ) and y=f^ ° ) (F.3) 

y o p 2 i N+I/ j y o pii N+u j 

are source terms for the pion condensate. This model can be extended to Nj pairs of 
quarks and conjugate quarks. Introducing P = iA + fiB and Q = iA^ + /IB' , we obtain 

Z™ Nf=2 (fi,p 1 ,p 2 ) = J dPdQ exp L 1 ^^} [ tr ( P tp + Qt Q) + tr(P Q + Qt P t)j I 

x p^det(ptp + p2 lAr+ ^ det (QtQ + ^ l7V )_ (F 4) 

To study the singular values of D{jl) we need to express the partition function in terms of 
the eigenvalues of P^P and Q^Q. Now, we notice that (F.4) is identical to eqs. (2.5)-(2.8) 
of [157]. 78 The representation of (F.4) in terms of the eigenvalues of P^ P and Q^Q is 
given in eq. (2.14) of that reference, and the microscopic spectral correlations were also 
computed in [157]. Hence, the microscopic correlations of the singular values of D(p) can 
be obtained immediately from the results of [157]. The mapping of RMT parameters to 
physical quantities can be worked out similarly to the /3 = 1 case in section 7. 

G Derivation of (7.52) 

From (7.51) we have for v = 

E'^ =2(y 4 $L) 2 |-ln/(z) with /(*)=/ dne z ^ +n "K (G.l) 
Performing a Hubbard-Stratonovich transformation gives 

f(z) = - dpdqe- p2+q2 dn e*V*<F»+v*) _ ( G2 ) 

J-oo JS 5 

78 It is interesting that the model of [157] is useful here since it applies to QCD with imaginary chemical 
potential, which is unrelated to QCD with real isospin chemical potential. 
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We now define a = (0, 2y/zp, 0, 2^/zq, 0, 0) and then rotate the coordinate system such that 
a is parallel to hi to obtain 



f dn e 2 ^ pn2+qn ^ = I 
Js$ Js 5 

/ d^sinV e acosv = ^rh{a) ■ (G.3) 
Jo a 



dn e ani with a = \a 

87T 2 f W , . „„„„,„ 87T 3 

~3 



After introducing polar coordinates for p and q we end up with 

4 3 /-oo -i? 2 

/(*) = / ^ -=- / 2 (2>/ii2) = vr 3 M(l, 3, z) , (G.4) 

z Jo u 

where M is Kummer's function (a.k.a. the confluent hypergeometric function). Using the 
recurrence relations for M we obtain 

M(l,3,z) = U e ^-l\, (G.5) 
and after performing the differentiation according to (G.l) we obtain (7.52). 
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